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ABSTRACT 


High  speed  data  conversion  is  required  to  conduct  digital  signal 
processing  on  wideband  analog  waveforms.  This  thesis  explores  the  performance 
characteristics  of  a  high  speed  optical  analog  to  digital  conversion  technique. 
This  technique  extends  the  resolution  of  an  integrated  optical  multi-interferometer 
analog  to  digital  converter  (ADC)  by  using  a  symmetrical  number  system  (SNS) 
encoding.  The  optical  SNS  ADC  can  directly  digitize  wideband  signals  with  a 
minimum  amount  of  circuit  complexity.  Using  an  8  bit  SNS  ADC,  the  results  are 
significantly  improved  with  a  modification  which  incorporates  a  parity  circuit,  and 
a  pulse  timing  circuit.  The  parity  circuit  decimates  the  possible  encoding  errors 
that  can  result  when  the  samples  lie  at  one  of  the  code  transition  points.  The 
results  of  this  decimation  process,  and  the  timing  of  the  corresponding  signals  are 
also  discussed.  Other  areas  such  as  the  integration  of  the  timing  processor  and 
folding  circuit,  are  identified  as  those  in  which  further  design  is  required. 
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L  INTRODUCTION 


A.  OVERVIEW  OF  ANALOG-TO-DIGITAL  CONVERTER 
L  Classincation  of  Anal(^  to  Digital  Converters 

Analog  to  digital  conversion  is  a  necessary  process  in  many  digital  electronic 
circuits.  This  process  converts  an  analog  signal  to  an  accurate  digital  number  proportional 
to  the  input  amplitude.  We  will  refer  to  devices  converting  analog  signals  to  digital 
representations  as  analog-to-digital  converters  (ADCs),  and  those  which  perform  the 
reverse  transformation  (digital  inputs  to  proportional  analog  signals)  as  digital-to-analog 
converters  (DACs).  Pulse-code  modulation  (PCM),  measurement  instmments  (multimeter, 
digital  memory  oscilloscopes  for  example),  signal-generation  and  processing  instruments 
(waveform  synthesizer  for  example)  are  all  common  applications  that  use  analog  to  digital 
conversion.[Ref.  1:  pp.  612-614] 

There  are  four  main  types  of  ADCs;  counter  (or  ramp)  ADC,  successive 
approximation  ADC,  flash  (or  parallel)  ADC,  and  integrating  ADC.[Ref.  2:  pp.  550-561] 
Each  type  has  its  advantages  and  disadvantages,  and  we  choose  one  of  these  types  based 
on  our  application.  The  counter  ADC  is  composed  of  four  parts:  the  comparator,  timing 
clock,  digital  counter,  and  DAC.  Figure  1  shows  a  block  diagram  of  a  counter  ADC.  We 
assume  that  the  initial  setting  of  the  counter  is  zero.  When  an  analog  sample  is  input  to  the 
comparator,  the  counting  process  is  initiated.  The  comparator  output  goes  into  the  digital 
counter  with  the  clock  pulse  gate.  Each  succeeding  clock  pulse  then  increments  the 
counter’s  binary  sequence  output.  The  DAC  translates  this  binary  code  to  an  analog  signal 
level.  The  comparator  then  compares  the  DAC  output  to  the  input  analog  sample.  When 
the  DAC  output  exceeds  the  level  of  the  analog  sample,  the  comparator  output  changes 
states,  stopping  the  digital  counter.  The  final  counter  state  is  the  desired  digital 
representation  of  the  analog  sample.  The  counter  ADC  is  the  simplest  type  to  understand. 
However,  it  is  relatively  slow  compared  to  the  others  because  of  the  many  comparisons 
which  must  be  performed  at  each  step.  A  counter  ADC  with  an  n-bit  output  may  require  as 

many  as  2“  - 1  comparisons  for  each  input  sample. 
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Digital 

output 


Figure  1.  Block  Diagram  of  a  Counter  ADC  from  Ref.  [2]. 
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Successive  approximation  is  another  common  general  purpose  ADC  type.  Figure  2 
shows  the  block  diagram  of  a  successive  approximation  ADC.  This  type  uses  control  logic 
and  an  output  register,  referred  to  as  a  successive  approximation  register  (SAR),  in  place  of 
the  clock  pulse  and  counter  of  the  counter  ADC.  The  comparator  compares  the  output  of 
the  DAC  to  the  input  analog  sample.  The  completion  of  conversion  is  triggered  by  a 
change  in  the  state  of  the  comparator.  The  diagram  of  both  the  counter  ADC  and  the 
successive  approximation  ADC  is  similar,  however,  the  logic  between  both  ADCs  is 
different. 

The  flow  chart  for  a  successive  approximation  ADC  is  shown  in  Fig.  3.  While  the 
counter  ADC  compares  the  next  state  of  converted-analog  signal  with  the  reference  signal, 
the  successive  approximation  ADC  compares  the  bit  information  from  the  most  significant 
bit  (MSB)  to  the  least  significant  bit  (LSB).  The  successive  approximation  ADC  thus 
requires  at  most  n  comparisons,  so  the  conversion  time  is  clearly  much  shorter.  A  typical 
conversion  time  of  the  successive  approximation  ADC  ranges  from  1  psec  to  50  p-sec. 

The  flash  ADC  has  the  shortest  conversion  time  of  any  available  type  of  ADC.  The 

block  diagram  of  a  flash  ADC  is  shown  in  Fig.  4.  An  n-bit  flash  ADC  requires  2"  - 1 

comparators.  The  input  analog  sample  is  applied  simultaneously  to  all  comparators.  The 
reference  voltage  for  each  comparator  is  set  through  a  common  resistor  ladder.  The 
decoder  converts  the  comparator  output  pattern  to  a  digital  word  representing  the  analog 
sample.  For  example,  if  we  have  seven  comparators  in  the  ladder  and  three  of  the  lower 
comparator  outputs  are  set  to  the  “ON”  state,  then  the  ladder  creates  a  digital  word 
“0000111.”  This  digital  word  representation  is  often  called  a  thermometer  code.  The 
decoder  converts  this  code  to  a  common  binary  representation  “Oil.”  For  a  decoder 

circuit,  a  read  only  memoty  (ROM)  is  widely  used.  Because  of  using  2“  -  1  comparators, 

we  do  not- need  to  set  up  a  feedback  loop.  We  need  only  one  comparison  of  the  analog 
sample  to  get  the  desired  digital  representation.  The  flash  ADC  is,  needless  to  say,  the 
fastest,  but  also  the  most  expensive  ADC.  In  addition,  the  number  of  bits  and  the  practical 
resolution  are  limited  because  of  the  many  comparators  loaded  in  parallel. 

The  integrating  ADC  is  used  for  q)ecial  applications.  This  type  uses  input  signal 
integration,  which  may  improve  the  signal-to-noise  ratio  (SNR)  for  certain  types  of  signals 
in  the  circuit.[Ref.  2:  pp.  557-561]  Thus,  using  an  integrating  ADC  may  help  suppress  the 
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Figure  2.  Block  Diagram  of  a  Successive  Approximation  ADC  from  Ref.  [2]. 
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Figure  3.  Flow  Chart  for  a  Successive  Approximation  ADC  from  Ref.  [2]. 
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F^ure  4.  Block  Diagram  of  a  Flash  (Parallel)  ADC  from  Ref.  [2]. 
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input  noise,  but  the  circuit’s  frequency  response  curve  distorts  the  spectrum  of  the  signal. 
Therefore,  the  integrating  ADC  is  used  only  for  dc  and  low  frequency  (slowly  varying) 
signals  to  preserve  the  integrity  of  the  signal.  Also,  this  type  of  ADC  is  relatively  slow  due 
to  the  precision  integrator  circuits.  For  example,  viith  the  Delta-Sigma  ADC  (A-E  ADC) 
which  uses  the  integrating  ADC  scheme,  the  conversion  rates  are  up  to  300  msec.  [Ref.  3] 
The  block  diagram  of  a  A-E  ADC  is  shown  in  Fig.  5.  The  input  voltage  of  the  A-E  ADC 
drives  the  integrator  and  is  compared  to  a  fixed  voltage  range  such  as  ground.  Depending 
on  the  comparator  output,  fixed  elements  of  charge  are  switched  into  the  summing  junction. 
A  counter  keeps  track  of  number  of  charge  pulses  switched  into  the  summing  junction,  and 
the  output  is  a  digital  word. 

2.  Current  Trend  of  Developing  ADCs 

ADCs  are  widely  used  to  translate  an  analog  signal  input  into  a  digital  signal  output. 
In  fields  such  as  digital  computing,  control  systems,  and  information  processing,  ADCs 
play  a  key  role.  Recent  advancements  in  scientific  technology  have  provided  for  a  shorter 
timing  response  in  signal  processing  therefore,  faster  sampling  rates  by  the  ADC  are 
required  fw  these  high  speed  applications.  As  mentioned  in  the  previous  section,  the  flash 
ADC  has  the  fastest  response  of  any  of  the  current  ADCs.  Some  research  papers  indicate 
an  effort  to  improve  the  sampling  rate  and  resolution  using  electronic  folding  devices.[Ref. 
4-6]  In  these  converters,  the  input  signal  is  folded  in  parallel  in  an  attempt  to  reduce  the 
complexity  of  the  amplitude  analyzing  function.  One  of  the  newest  high  speed  folding 
ADC  architectures  is  the  integrated  optical  guided  wave  (lOGW)  technique.  [Ref.  7-13] 
The  lOGW  technique  can  handle  higher  sampling  rates  than  conventional  ADCs  and  can  be 
used  for  high  speed  applications. 


B.  lOTEGRATED  OPTICAL  SNS  ADCS 
L  Symmetrical  Number  System 

The  symmetrical  number  system  (SNS)  is  composed  of  a  number  of  pairwise 
relatively  prime  (PRP)  moduli  mi.  The  integers  within  each  SNS  modulus  are 
representative  of  a  symmetrically  folded  waveform  with  the  period  of  the  waveform  equal 
to  twice  the  PRP  modulus.  For  m  given,  the  integer  values  within  twice  the  individual 
modulus  Xm  are  given  by: 
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(1) 


The  SNS  folding  waveform  is  symmetrical  about  the  midpoint  and  is  compatible  with  other 
folding  waveforms.  Due  to  the  presence  of  ambiguities,  the  integers  within  Eq.  1  do  not 
form  a  complete  system  of  length  2m  by  themselves.  It  is  well  known  however,  that  the 
inclusion  of  additional  redundant  moduli  can  effectively  detect  and  correct  errors  within  a 
residue  number  system  (RNS)  representation  of  a  number.  The  SNS  formulation  is  based 
on  a  similar  concept  which  allows  the  ambiguities  to  occur.  The  ambiguities  that  arise 
within  the  SNS  are  resolved  by  using  various  arrangements  of  the  SNS  moduli.  By 
considering  the  derived  moduli  arrangements,  the  SNS  is  rendered  a  complete  system 
having  a  one-to-one  correspondence  with  the  RNS.  For  N  equal  to  the  number  of  PRP 
moduli,  the  dynamic  range  M  of  the  system  is  as  follows  [Ref.  1 1]: 

N 

M=nw..  (2) 

i=l 

where  M  is  also  the  position  of  the  first  repetitive  moduli  vector. 

2.  Integrated  Optical  SNS  ADC 

The  most  efficient  (sampling)  ADCs  demonstrated  to  date  use  lOGW  technology  (in 
lithium  niobate  (LiNbOs),  gallium  arsenide  (GaAs),  or  indium  phosphide  (InP)). 
Symmetrical  folding  of  the  input  signal  is  realized  with  Mach-Zehnder  Intaferometer  (MZI) 
waveguide  modulators  arranged  in  a  parallel  configuration.  [Ref.  7-10]  These  modulators 
use  the  linear  Pockels  effect,  and  provide  a  convenient  method  for  coupling  a  wideband 
electrical  signal  into  an  optical  processing  system  through  the  modulator  electrodes.  This 
architecture,  often  called  a  Taylor  scheme,  makes  use  of  the  periodic  dependence  of  the 
modulator  optical  output  on  the  applied  voltage  to  the  electrodes  and  the  electrode  lengths. 
The  ultimate  bandwidth  of  the  signal  that  can  be  digitized  is  determined  by  the  sampling 
pulse  width  and  rate,  the  modulator  bandwidth,  the  detector  bandwidth,  and  the  speed  of 
the  output  digital  circuitry.  One  of  the  major  limitations  for  the  Taylor  scheme,  however,  is 
the  required  doubling  of  the  electrode  lengths  for  the  increasing  number  of  bits.  For 
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example,  an  n-bit  ADC  using  a  Taylor  scheme  would  require  electrode  lengths  of  /,  21, 41, 
2”"*  •/  where  the  length  of  the  longest  modulator  is  2""*  •/.  Here  /  stands  for  the 

shortest  electrode  length.  This  electrode-length  doubling  requirement  results  in  a  high 
device  capacitance,  and  ultimately  constrains  the  achieved  resolution  to  about  2  to  4  bits, 
which  is  not  sufficient  for  many  applications.[Ref.  10] 

Recently,  a  technique  that  extends  the  resolution  of  an  integrated  optical 
multi-interferometer  analog-to-digital  converter  was  described.  [Ref.  11-13]  The  optical 
output  waveform  for  each  interferometer  is  symmetrically  folded  at  twice  a  proper 
modulus.  A  small  comparator  ladder  mid-level  quantizes  each  interferometer's  detected 
ouq)ut  to  encode  the  analog  signal  in  a  SNS  format.  By  incorporating  the  SNS  encoding, 
resolution  greater  than  1-bit  per  interferometer  can  be  provided. 

By  applying  Eq.  1  and  Eq.  2,  the  SNS  can  serve  as  a  source  for  resolution 
enhancement  in  an  ADC  by  decomposing  the  analog  amplitude  analyzing  function  into  a 
number  of  parallel  sub-operations  (moduli)  that  are  of  smaller  computational  complexity. 
Each  sub-operation  for  a  different  modulus  requires  only  a  precision  based  on  that 
modulus.  A  much  higher  resolution  is  achieved  after  the  results  of  these  low  precision 
sub-operations  are  recombined.  That  is,  the  resolution  of  each  interferometer  is  greater 
than  1-bit  per  interferaneter.  An  original  schematic  diagram  of  a  SNS  ADC  architecture  is 
shown  in  Fig.  6.  The  input  signal  is  folded  in  parallel  with  each  folding  period  equal  to 
twice  a  particular  modulus.  The  folded  waveform  at  the  output  of  each  folding  circuit  is 
mid-level  quantized  with  a  small  comparator  ladder  to  encode  the  input  signal  into  the  SNS 
format  An  encoder  then  converts  the  SNS  representation  to  a  more  familiar  digital  ouq)ut 
such  as  a  binary  representation.  With  the  SNS  encoding  any  combination  of  folding 
periods  and  comparator  arrangements,  any  analysis  can  be  exact 

The  optical  output  power  from  a  single  guided- wave  interferometer  is  symmetrical 
and  periodic,  and  can  thus  be  used  to  implement  each  folding  circuit  (modulus)  in  the  SNS 
ADC.  The  interferometer  normalized  output  is  a  function  of  both  the  sampled  analog 
voltage  V  and  the  modulus  m,-  as 


/(v..)  =  sin  (-^)  ie{l,2,...,N} 


(3) 
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Figure  6.  Original  System  Diagram  of  an  SNS  ADC 
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where  v,-  is  the  modulus  dependent  attenuated  input  signal,  and  N  is  the  number  of 

modulators.  The  attenuation  factor  depends  on  the  modulus  in  such  a  way  that  identical 
interferometers  can  be  used.  Specifically,  we  can  take 


m. 


v.= 


■V 


(4) 


where  V  is  the  input  signal,  and  V„  is  the  half  wave  voltage  for  the  MZI.[Ref.  14:  pp.  324- 
330]  The  rrimin  is  the  minimum  modulus  of  the  set  needed  to  satisfy  Eq.  4. 

The  folding  circuit  for  each  modulus  requires  m,-  - 1  comparators  at  the  ouqrut  of  the 

detector.  With  the  waveform  given  by  Eq.  3,  the  normalized  threshold  values  for  the 
comparators  within  each  modulus  m,-  are 


where 


T  - 

V 


«  ^  V  • 

sin^y-^) 


(5) 


v^.= 


j  6{l,2,...,m.-l}. 


(6) 


That  is,  the  comparator  thresholds  are  tailored  to  the  interferometer  output  waveform.  The 
size  of  the  LSB  for  the  SNS  ADC  is  obtained  from  Eq.  4  and  Eq.  6,  under  the  condition 
that  nti  =  rrimin. 


V  = 

LSB 


(7) 
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and  therefore  a  larger  minimum  modulus  leads  to  higher  resolution. 

Table  1  shows  several  possible  SNS  ADC  configurations.  Systems  with  resolution 
on  the  order  of  seven  to  eleven  bits  are  shown.  Detailed  are  the  moduli  corresponding  to 
each  interferometer  with  the  dynamic  range,  the  total  number  of  comparators  required,  and 
the  maximum  number  of  comparators  loaded  in  parallel  at  a  detector  output  (fan  out). 

Another  device  constraint  is  the  maximum  voltage  that  may  be  applied  to  the 
electrodes.  Applied  bipolar  voltage  beyond  the  rated  maximum  (minimum)  will  spark 
across  the  electrode  structure  and  damage  the  device.  The  specifications  of  and  Vmax 
reveal  the  maximum  number  of  folds  available  from  the  device.  An  input  voltage  which 
makes  one  complete  fold  is  2  .  The  maximum  number  of  folds  F,fuix  available  from  the 

device  is  therefore 


Pmax  = 


2v. 


2V 


(8) 


The  smallest  modulus  within  the  SNS  system  requires  the  largest  number  of  folds 
to  instrument  the  dynamic  range.  Note  that  each  fold  in  the  minimum  modulus  covers 
2/WOTin  SNS  ADC  states.  The  largest  number  of  folds  required  from  an  interferometer  in  a 
5-bit  SNS  ADC  is  therefore 


F 


req 


(9) 


C  PRINCIPAL  CONTRIBUTIONS 

In  this  thesis,  MATLAB™  programming  code  of  the  SNS  ADC  is  provided  as  a 
simulation  of  the  theoretical  system  behavior.  The  initial  program  offers  the  normalized 
threshold  level  and  the  resistance  value  information  in  the  ladder  circuit  of  the  SNS  ADC, 
and  the  digital  word  conversion  between  the  thermometer  code  and  the  binary  code  for  an 
arbitrary  modulus.  Furthermore  an  8-bit  SNS  ADC  system  simulation  with  a  parity  circuit 
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Interferometer 

Dynamic  Range 

Total  Number 

Max  Number 

Number 

ModuU 

of  Comparators 

Lxmded  in  Parallel 

of  Bits 

M 

L 

K 

B 

3,4,5 

60 

9 

4 

5 

11,12 

132 

21 

11 

7 

5,7,8 

280 

17 

7 

8 

9,10,11 

990 

27 

10 

9 

13,14,15 

2730 

39 

14 

11 

Table  1.  Possible  SNS  ADC  Configurations 
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to  discard  encoding  errors  is  described.  MATLAB  is  also  used  in  this  simulation.  The 
program  computes  the  normalized  threshold  levels  of  the  parity  circuit,  the  resistance 
values,  and  each  resistance  value  in  the  ladder  circuit.  The  parity  circuit  reduces  the 
encoding  error  problems  of  the  original  SNS  ADC,  and  this  is  measured  using  the  logic 
analyzer.  Comparison  and  discussion  of  the  results  between  the  actual  SNS  ADC  output 
using  Lab  VIEW  generated  inputs  and  the  output  computed  by  the  MATLAB  simulation  is 
accomplished.  Decimation  using  the  parity  circuit  are  also  described. 

Because  of  multi-interferometer  operation,  timing  is  needed  to  maintain  proper 
operation  of  the  SNS  ADC.  Therefore  a  timing  processor  is  developed,  and  the  circuit  is 
constructed.  Each  component  of  the  timing  processor,  and  the  timing  characteristics  of  the 
processor  are  described. 

D  THESIS  OUTLINE 

Chapter  n  of  this  thesis  consists  of  an  overview  of  MATLAB  followed  by  a 
description  of  the  initial  simulation  development.  Next,  the  extended  development  of  the 
MATLAB  program  to  include  the  parity  analysis  is  discussed.  Finally,  total  system 
simulation  is  obtained  with  MATLAB  and  documented.  Chapter  m  provides  the  details  on 
the  timing  processor  of  the  system,  and  a  description  of  the  components  considered. 
Chapter  IV  shows  some  experimental  results  of  the  actual  circuit,  and  discusses  the 
difference  between  the  results  simulated  by  MATLAB  and  the  results  of  the  actual  circuits. 
Following  the  summary  in  Chapter  V,  the  Appendix  provides  the  MATLAB  code  that 
computes  the  design  and  simulates  the  whole  SNS  ADC  system  input  This  input  uses  the 
same  signal  generating  methods  as  Lab  VIEW.  Each  program  in  the  Appendix  is  discussed 
in  Chapter  H. 
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n.  SIMULATION  OF  THE  SNS  ADC  USING  MATLAB™ 


A.  OVERVIEW 

In  optical  integrated  circuits,  the  linear  electrooptic  effect  (or  the  Pockels  effect)  is 
widely  used,  especially  for  light  modulation  and  switching.  An  index  change  is  induced  in 
the  crystals  via  the  electrooptic  effect  when  an  electric  field  is  applied  to  the  crystal.  When 
the  index  varies  linearly  with  the  applied  electric  field,  it  is  called  the  linear  electrooptic 
effect.  Since  the  linear  electrooptic  effect  is  found  in  crystals  without  an  inversion 
symmetry,  such  as  LiNbOs,  the  sign  of  the  induced  index  change  depends  on  the  polarity 
of  the  voltage  applied  to  the  crystal.  The  MZI  is  one  of  the  popular  optical  computing 
devices  that  use  the  linear  electrooptic  effect  to  induce  phase  modulation  in  the  optic 
signals.  Equation  3  shows  the  output  of  the  MZI  which  is  a  sinusoidal  function.  Each 
MZI  folds  the  input  signal  before  being  digitized  by  high  speed  comparators. 

The  electronic  part  of  the  original  SNS  ADC  is  shown  in  Fig.  7.  It  includes  a 
comparator  ladder  circuit,  priority  encoder,  and  EPROM  (erasable  programmable  ROM). 
The  EPROM  is  similar  to  the  programmable  logic  array  (PLA).  The  EPROM  is  a 
programmable  and  erasable  device,  which  converts  the  thermometer  code  into  a  binary 
representation.  In  an  ADC  SNS  design,  the  threshold  voltage  levels  of  the  comparator 
ladder  for  each  modulus  need  to  be  derived.  However,  the  threshold  voltage  level  for  each 
modulus  is  different.  One  of  the  reasons  for  the  different  voltage  levels  is  the  non-linear 
characteristic  of  the  modulator,  and  the  differences  in  the  folding  periods.  Also  the 
threshold  levels  depend  on  each  modulus  in  the  SNS.  According  to  the  SNS  sonicture,  m- 
1  threshold  values  are  needed  for  mod  m  and  these  threshold  voltages  can  be  expressed  by 
applying  Eq.  3.  The  characteristic  curve  of  the  modulator  of  each  MZI  is  the  same,  and  we 
can  get  the  value  of  the  predicted  threshold  voltage  for  each  modulus  using  computer 
simulation. 

MATLAB  is  a  technical  computing  environment  for  high-performance  numeric 
computation.  It  integrates  numerical  analysis  and  matrix  operation,  and  its  interface  is  very 
easy  to  use.  MATLAB  can  also  be  configured  for  various  computer  platforms.  The 
MATLAB  program  developed  in  this  thesis  is  discussed  in  the  following  section. 
MATLAB  is  a  trademark  of  The  Math  Works,  Inc. 
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Figure  7.  The  Electronic  Part  of  the  SNS  ADC 
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B.  DEVELOPMENT  OF  MATLAB  PROGRAM 


L  Computing  the  Threshold  Values  for  Each  Modulus 

The  program  to  calculate  the  normalized  threshold  voltages  (MATLAB  code)  is 
shown  in  Appendix  A.  The  normalized  value  of  the  threshold  voltage  for  each  comparator 
in  mod  m  can  be  obtained  using  the  program  called  Vth_calc.m.  This  program  applies  Eq. 
3,  which  details  the  relationship  between  the  intensity  of  the  MZI  output  and  its 
corresponding  input  voltage.  The  m-1  threshold  values  in  mod  m  are  calculated,  and  each 
threshold  value  depends  on  the  input  voltage  to  the  SNS  ADC  (D).  An  example  drawing 
showing  the  threshold  voltages  (mod  9)  is  shown  in  Fig.  8.  In  this  case,  eight  threshold 
values  are  needed  and  each  value  can  be  calculated  using  the  first  half  period  of  the  MZI 
ouq)ut.  The  small  circles  on  the  graph  show  the  threshold  voltage  levels  for  each  input 
voltage  index.  The  x-axis  shows  the  index  number  for  the  input  voltage  D,  and  the  y-axis 
shows  the  normalized  intensity  of  the  MZI  output.  Actually,  we  cannot  measure  the  light 
intensity  directly,  but  we  can  measure  the  voltage  as  an  intensity  related  value  using  a  photo 
detector.  If  we  know  the  maximum  voltage  of  the  MZI  output  max(VMZ^oMr).  we  can  get 
each  threshold  voltage  as; 


Vif/i„y=max(Vj 


MZJl  out 


)T..  ie{l,2,...,iV} 

y  6{l,2,...,m-l} 


(10) 


The  Vth_calc.m  program  can  also  be  modified  for  the  general  case,  in  which  the  minimiim 
intensity  of  MZI  output  is  not  equal  to  zero.  In  this  case,  the  threshold  voltage  can  be 
calculatetTusing  the  minimum  and  maximum  voltage  of  the  photo  detector  output; 

(max(V«z,,^)-min(Vj»^,^))  T.  . 

ie{l,2,...,iV}  (11) 

y  6  {1,2,...,  m-1} 

The  SNS  ADC  uses  multiple  interferometers,  and  therefore  calculating  the  threshold 
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Figure  8.  Example  of  Normalized  Threshold  Values  vs.  Input  Voltage  Index  for 
mod  9 
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voltages  for  the  various  moduli  is  needed  for  each  system.  The  Vth_calc.m  handles  this 
process  also.  This  program  also  returns  the  difference  of  neighboring  normalized  threshold 
values  (dVth).  Vth  and  dVth  values  are  required  for  building  up  the  comparator  ladder  of 
the  flash  ADC  scheme  for  each  modulus. 

Appendix  B  shows  the  program  named  thmtobin.m  which  converts  the 
thermometer  code  generated  by  the  priority  encoder  into  the  binary  code.  This  program  is 
needed  for  confirming  operation  of  the  electronic  part  of  the  SNS  ADC.  The  user  enters 
the  information  for  each  modulus,  and  the  computer  returns  the  possible  states  of  the 
comparator,  thermometer  code  representation,  and  the  binary  information.  Using  these 
two  programs  Vth_calc.m  and  thmtobin.m,  we  can  get  the  information  needed  to  constract 
the  SNS  ADC  ciicuit.[Ref.  15] 

2.  Creating  the  Folding  Waveform  Inputs  Using  LabVIEW™ 

Lab  VIEW™  is  a  program  development  application  which  uses  a  graphical 
programming  language.  Virtual  instmments  (Vis)  are  one  of  the  LabVIEW’s  powerful 
features.  Using  the  Vis,  we  can  develop  a  specific  program  for  data  acquisition  and 
instrumental  control.  LabVIEW  is  a  trademark  of  National  Instramenls  Corporation.  In 
the  early  stage  of  the  SNS  ADC  circuit  development,  LabVIEW  was  used  to  simulate  the 
input  signal  to  the  comparators.fTRef.  16,17]  For  the  folding  circuits,  there  are  three  design 


rules,  [Ref.  11] 

Rule  (1) 

M  <2  X  (m  X  (number  of  folds) 

(12) 

Rule  (2) 

2®-i 

i 

(13) 

Rule  (3) 

V 

T?  ^ 

^  rtq  -  y 

(14) 

where  B  is  the  number  of  bits,  Freq  is  the  required  number  of  the  folds,  Vmax  is  the 
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maximum  input  voltage,  Vjc  is  a  half  wave  voltage  for  the  MZI,  and  M  is  the  dynamic 

range  of  the  SNS  ADC.  Using  these  relationships,  the  folding  simulation  program  of 
LabVIEW  is  set  up.[Ref.  16]  The  goal  of  the  initial  effort  is  to  simulate,  and  compare  the 
results  generated  using  LabVIEW  and  the  digital  hardware  to  the  output  of  the  MATLAB 
simulation  program.  Figure  9  shows  the  folding  waveform  ouq)ut  of  an  8-bit  SNS  ADC. 
In  this  case  mi  =  9,  m2  =  10,  and  m3  =  11.  This  figure  was  obtained  by  applying  Eq.  3 
The  x-axis  shows  the  input  voltage  index  and  the  y-axis  shows  the  normalized  value  of  the 
MZI  output.  Note  that  the  initial  state  of  each  MZI  output  is  set  to  the  minimum  value  of 
the  folding  waveform. 

In  LabVIEW,  files  called  JEFFSTEP.VI  and  CRAIG7.VI  are  used  to  simulate  the 
folding  circuit  output.[Ref.  16,17]  LabVIEW  creates  the  step  generation  of  a  sine  wave 
with  the  program  files  above.  In  the  8-bit  SNS  ADC  design,  we  know  that  the  MZI  which 
is  used  in  the  actual  ADC  SNS  circuit,  has  a  folding  period  of  2Vn  or  4.507  V  in  mod  9, 

and  32  V  is  supplied  as  v^ax-  This  maximum  voltage  gives  a  sufficient  number  of  folding 

waves.  The  number  of  folds  for  mod  9  is  then  calculated  using  Eq.  14,  and  the  number  of 
folds  for  the  other  moduli  are  calculated  by  using  the  ratio  of  the  lowest  modulus  to  the 
modulus  of  concern.  In  this  case  9/10  and  9/1 1  are  used  respectively  for  mod  10  and  mod 
1 1.  Let  Fm  be  the  number  of  folds  for  mod  m.  Then  F9,  F  jo,  and  Fn  are  calculated  firom 
Eq.  8  as  follows: 


''10  =  12.78  (16) 

P.i  (17) 

Using  the  cycle  ratio  for  each  folded  output  waveform,  the  frequency  ratios  are  calculated. 
The  files,  JEFFSTEP.VI  and  CRAIG7.VI  use  110,  99,  and  90  cycles  for  mod  9,  10  and 
11  respectively  in  the  8-bit  design.  The  number  of  points  of  observation  (#obs)  is 
determined  by  following[Ref.  16]: 
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Figure  9.  The  Folding  Waveform  Output  of  an  8-bit  SNS  ADC 
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5,000  is  initially  chosen  as  the  number  of  generation  points  on  both  JEFFSTEP.VI  and 
CRAIG7.VI  to  represent  the  whole  dynamic  range.  Since  the  number  of  cycles  associated 
with  the  lowest  modulus  (mod  9)  is  1 10,  5,(X)0/1 10  corresponds  to  one  fold  of  mod  9. 
Therefore,  the  number  of  points  of  observation  for  the  8-bit  design,  which  requires  only 
256  states,  is  646  by  Eq.  15,  and  Eq.  18  as  follows: 


Number  of  generation  points  ^  ^  5,0(X) 

lowest  modulus  cycles  ®  ~  110 

=  645.5 


(19) 


The  step  input  increment  (size  of  the  least  significant  bit)  V  lsb  is  determined  to  be 
0.250  V  by  Eq.  7.  The  size  of  the  least  significant  bit  for  each  modulus  is  the  same.  In 
order  to  compare  with  the  LabVIEW  results,  we  call  the  number  of  observation  points  per 
input  voltage  index  the  number  of  sample  points.  The  number  of  sample  points  (#sampie) 
must  be  calculated  and  should  also  be  the  same  for  each  modulus,  since  each  modulus  has 
same  sample  points  in  a  system.  This  number  of  sample  points  can  be  computed  from 


f#  )= 
'  sample' 


number  of  generation  points 
number  of  cycle  x  2  x  m 


For  example  using  5,0(X)  points,  110  cycles,  and  m  =  9,  Eq.  20  provides  #sampie  =  2.5. 

The  MATLAB  code  is  programmed  using  the  information  above.  Appendix  C 
shows  the  simulating  LabVIEW  program  file.  This  MATLAB  file  named  Lab_simu.m  uses 
the  following  input  information. 

•  Number  of  generating  points 

•  Number  of  points  of  observation 

•  Number  of  cycles  of  each  modulus 

•  Number  of  each  modulus 
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This  infomiation  provides  die  folding  circuit  information  as 

•  The  number  of  folds  of  mod  m,  Fm 

♦  The  number  of  sanqile  points  between  two  thresholds  of  each  comparator 

*  The  tiiermometer  code  for  each  sample  point  of  each  modulus 

•  Predicted  normalized  LabVIEW  output  at  each  sample  point 

Input  information  for  this  program  should  be  prepared  in  advance.  For  example,  an  8-bit 


SNS  ADC  design  using  mod  9,  10,  and  11  with  5,000 
information  is  as  follows: 

generating  points,  the  input 

• 

Number  of  generating  points  : 

5,000 

• 

Number  of  points  of  observation  : 

647  (from  0  to  646) 

• 

Number  of  cycles  of  each  modulus  : 

110,99,90 

• 

Number  of  each  modulus  : 

9,  10,  11 

3. 

Simulation  Results 

After  programming  the  file  Lab_simu.m,  the  output  of  both  the  MATLAB 
simulation  and  LabVIEW  are  compared.  The  LabVIEW  program  named  CRAIG7.VI 
creates  the  output  value,  and  eidiCT  displays  or  prints  out  the  numbers,  this  file  uses  5.0  V 
as  a  maximum  voltage  of  the  folding  circuit,  so  each  voltage  value  of  the  index  obtained  by 
the  MATLAB  simulation  must  be  multiplied  by  five.  The  output  of  the  MATLAB 
simulation  always  provides  a  normalized  value.  CRAIG7.VI  calculates  three  decimal 
places  instead  of  the  four  decimal  places  calculated  with  the  MATLAB  simulation.  Table  2 
shows  the  comparison  of  some  of  the  output  values  generated  with  CRAIG7.VI  and 
Lab_simu.m. 

As  a  result,  the  MATLAB  simulation  can  simulate  the  LabVIEW  generated  signals. 
This  program  can  also  simulate  the  thermometer  code  output  for  each  modulus.  These 
outputs  are  very  helpful  to  confirm  the  output  states  of  each  modulator.  On  an  actual  8-bit 
SNS  ADC  board,  three  LED  displays  are  used  for  confirming  these  ou^ut  states.  The 
MATLAB  simulation  can  also  be  used  to  create  an  arbitrary  configuration  of  an  SNS  ADC 
design.  The  file,  Lab_simu.m,  is  based  on  the  interferometer  combination  of  mod  9, 10, 
and  11  with  5,000  generating  points,  and  we  can  apply  this  file  to  a  new  combination 
sin:q)ly  by  changing  the  input  information. 
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Index 

LabView 
(mod  9) 

LabView 
(mod  10) 

LabView 
(mod  11) 

MATLAB 
simulaion 
(mod  9) 

MATLAB 
simulaion 
(mod  10) 

MATLAB 
simulaion 
(mod  11) 

1 

0 

0 

0 

0 

0 

0 

2 

0.024 

0.019 

0.016 

0.0238 

0.0193 

0.0160 

3 

0.095 

0.064 

0.0949 

0.0770 

4 

0.212 

0.143 

5 

0.373 

0.303 

0.3725 

0.3032 

0.2515 

6 

0.574 

0.468 

7 

0.812 

0.665 

0.554 

0.8117 

0.6647 

0.5538 

8 

1.082 

0.890 

0.743 

1.0818 

0.8896 

0.7434 

9 

1.379 

1.139 

0.955 

1.3790 

1.1393 

0.9554 

10 

1.698 

1.410 

1.187 

1.6976 

1.4100 

1.1871 

11 

2.032 

1.698 

1.436 

2.0315 

1.6976 

1.4356 

12 

2.374 

1.998 

1.698 

2.3744 

1.9976 

1.6976 

13 

2.720 

2.305 

1.970 

2.7196 

2.3054 

1.9700 

14 

3.061 

2.616 

2.249 

3.0607 

2.6162 

2.2491 

15 

3.391 

2.925 

2.531 

3.3910 

2.9252 

2.5314 

16 

3.704 

3.228 

2.813 

3.7044 

3.2276 

2.8133 

17 

3.091 

3.5188 

3.0912 

18 

4.257 

3.794 

3.362 

4.2566 

3.7942 

3.3616 

19 

4.485 

4.050 

3.621 

4.4850 

4.0496 

3.6210 

20 

4.675 

4.281 

3.866 

4.6755 

4.2810 

3.8660 

21 

4.824 

4.485 

4.094 

4.8244 

4.4850 

4.0936 

iOil 

4.929 

4.658 

4.301 

4.9291 

4.6582 

4.3008 

23 

4.987 

4.798 

4.485 

4.9874 

4.7981 

4.4850 

24 

4.998 

4.902 

4.644 

4.9982 

4.9025 

4.6438 

25 

■Sttil 

4.775 

4.9614 

4.9697 

4.7753 

26 

4.878 

4.999 

4.878 

4.8776 

4.9988 

4.8776 

Table  2.  Comparison  of  Program  Output  Value  with  LabVIEW  and  MATLAB 
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C  EXTENDED  DEVELOPMENT  OF  MATLAB  PROGRAM 
L  Overview  of  the  Parity  Circuit 

The  SNS  ADC  simulation  program,  developed  in  the  previous  section,  is  based  on 
the  original  SNS  ADC  design  shown  in  Fig.  6  in  Chapter  1.  This  original  design, 
however,  results  in  a  number  of  encoding  errors. 

In  the  original  design  of  the  SNS  ADC,  it  is  assumed  that  the  state  of  each  modulus 
changes  simultaneously.  Figure  10  shows  details  of  a  folding  waveform  of  an  8-bit  SNS 
ADC.  The  important  thing  about  the  SNS  ADC  is  that  the  folding  waveforms  for  each 
modulus  pass  through  the  threshold  values  simultaneously  across  all  channels.  This  is 
demonstrated  by  the  vertical  line  drawn  at  the  input  voltage  index  of  six.  Simultaneous 
transition  of  all  modes  from  one  state  to  another  is  the  ideal  situation.  However,  this  may 
not  happen  in  an  actual  circuit  because  of  the  inaccuracies  of  the  matching  voltages.  When 
the  folding  waveforms  do  not  pass  through  their  matching  threshold  levels  simultaneously, 
an  encoding  error  results. 

To  eliminate  the  possibility  of  these  encoding  errors,  a  parity  circuit  is  used  to  set 
up  a  band  around  each  code  transition  point  This  is  illustrated  in  Rg.  1 1  for  one  Areshold 
level.  The  circuit  modification  requires  a  few  additional  comparators  in  the  smallest 
channel.  That  is,  2m;  comparators  are  used  instead  of  m;  -1.  Of  these,  2m;  -  2  would  be 

paired,  one  just  below  the  proper  comparator  threshold  level  and  the  other  just  above.  Of 
the  remaining  two,  one  would  be  just  above  the  minimum  modulation  depth,  and  the  other 
threshold  just  below  the  maximum  detector  output.  A  schematic  diagram  of  the  SNS  ADC 
including  the  parity  circuit  is  shown  in  Fig.  12.  The  operation  of  this  parity  circuit  is  as 
follows: 

Assume  that  a  signal  arrives  and  all  comparators  below  a  certain  level  are  on.  This 
is  as  it  should  be.  For  each  incoming  signal,  if  it  turns  on  an  even  number  of  comparators, 
0,2,4,  ...,  m;,  the  signal  is  rejected.  If  it  turns  on  an  odd  number  of  comparators,  it  is 

accepted.  In  this  manner  we  place  a  narrow  band  around  each  voltage  level  that 
corresponds  to  a  switching  point  For  all  moduli,  the  comparator  levels  are  adjusted  so  that 
they  cross  their  respective  folding  waveforms  within  these  narrow  bands.  Any  input  signal 
within  these  narrow  rejection  bands  can  be  easily  discarded.  The  parity  circuit  operation 
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Figure  11.  The  Concept  of  Parity  Circuit 
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MZI  No.l 


D 


Figure  12,  Block  Diagram  of  the  SNS  ADC  Including  the  Parity  Circuit 
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region  of  mod  9  is  shown  in  Fig.  13.  The  gray  shaded  region  indicates  the  rejection  area. 
Samples  in  this  region  will  be  ignored.  The  unshaded  areas  indicate  the  acceptance  region. 
Legitimate  state  combinations  will  be  accepted  as  a  valid  ou^ut  The  dotted  horizontal  line 
is  the  original  threshold  level  of  mod  9.  The  solid  line  above  and  under  the  dotted  line 
shows  the  iqrper  and  lower  threshold  levels  for  the  parity  comparators. 

2.  Parity  Circuit  Threshold  Voltages 

To  set  up  the  parity  circuit  threshold  levels,  we  must  select  reasonable  rejection 
(decimation)  regions.  First  we  need  to  get  Vlsb-  Once  Vlsb  is  obtained,  the  width  of  the 
decimation  region  can  be  calculated  from  the  user’s  width  specified  as  a  percentage  of  the 
LSB  width. 


(decimation  width) = 


(%desired) 

ioo 


LSB 


(21) 


The  upper  and  lower  parity  circuit  input  voltage  for  input  voltage  index  k  are  chosen  as 
follows; 

(upper  parity);^  ic  ~  ^in,  k  y  (decimation  width)  (22) 

(lower  parity);^  k  ~  ^in.  k  “  ^  (decimation  width)  (23) 


where  is  an  input  voltage  corresponding  to  the  input  voltage  index  k.  Finally  we  can 
compute  the  2(m- 1 )  +2  parity  circuit  normalized  threshold  levels  as 


(upper  parity)  sin^  (upper  parity) 

(24) 

(lower  parity)  =  sin^  (lower  parity)  . 

(25) 

Figtne  14  shows  an  example  of  this  parity  circuit  computation.  In  this  example,  mod  3  is 
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NORMALIZED  VALUE  OF  MZI  OUTPUT 
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Figure  14.  Example  of  the  Parity  Circuit  Threshold  Computation  (mod  3) 
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used  for  the  smallest  channel,  and  the  decimation  width  is  60%  of  Vlsb-  The  larger  the 
percentage  of  decimation  width  chosen,  the  wider  the  threshold  levels  between  the  parity 
and  rejection  range. 

The  MATLAB  program,  idlprty.m,  covers  these  processes.  Idlprty.m  is  given  in 
Appendix  Section  D.  The  process  to  set  up  the  parity  circuit  threshold  levels  is  similar  to 
the  program  Vth_calc.m.  Desired  percentage  of  the  rejection  region  and  the  modulus 
number  of  the  least  significant  bit  are  needed.  Using  the  decimation  width  and  the  modulus 
corresponding  to  the  smallest  channel,  the  program  first  computes  the  input  voltage  indices 
and  Vlsb,  and  the  original  threshold  values.  The  threshold  levels  for  the  parity  circuit  are 
calculated  in  the  next  stage,  and  finally,  the  output  values  of  the  parity  circuit  are  computed. 
Figure  15  shows  the  threshold  levels  for  a  smallest  channel  of  mod  9  with  a  20% 
decimation  width. 

3.  Simulating  the  System  with  the  Parity  Circuit 

The  final  goal  of  the  MATLAB  SNS  ADC  simulation  is  to  simulate  the  complete 
system  using  the  LabVIEW  signals.  Required  output  includes: 

1 .  Threshold  level  for  each  modulus 

2.  The  input  voltage  samples  generated  by  LabVIEW 

3.  Ouq)ut  state  related  to  each  input  voltage  sample 

4.  Parity  circuit  threshold  levels  of  the  smallest  modulus 

5.  The  closest  states  and  the  difference  for  each  sample  output 

6.  Parity  bit  information  (accept  /  reject) 

The  first  four  items  are  completed  by  Vth_calc.m,  Lab_simu.m,  and  idlprty.m.  The 
MATLAB  file,  sys_simu.m  is  shown  in  Appendix  Section  E  and  runs  all  three  of  these 
files.  In  order  to  meet  the  fifth  requirement  above,  the  following  technique  is  used. 
Initially,  two  vectors  are  created.  The  vector  named  AdiffmVout  shows  the  difference 
between  the  ouq)ut  of  each  sample  index  and  the  closest  threshold  just  below  the  sample  for 
a  positive  slope  of  the  folding  waveform,  or  the  difference  between  the  output  of  each 
sample  index  and  the  closest  threshold  just  above  the  sample  for  a  negative  slope  of  the 
folding  waveform  of  mod  m.  BdiffmVout  shows  the  difference  between  each  sample 
index  and  its  closest  threshold  just  above  the  sample  for  a  positive  slope  folding  waveform, 
or  the  difference  between  each  sample  index  and  its  closest  threshold  just  below  the  sample 
for  a  negative  slope  folding  waveform  of  mod  m.  Once  obtained,  these  two  vectors  are 
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NORMALIZED  OUTPUT  VOLTAGE 


Normalized 
Parity  Threshold 
levels  for  mod  9 

0.0003 

0.0245 

0.0364 

0.1060 

0.1284 

0.2350 

0.2653 

0.3960 

0.4304 

0.5696 

0.6040 

0.7347 

0.7650 

0.8716 

0.8940 

0.9636 

0.9755 

0.9997 

(In  case  of  20  % 
decimation  width) 


Figure  15.  Example  of  the  Parity  Circuit  Threshold  Computation  (mod  9) 
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compared  with  each  other.  The  smallest  value  is  chosen,  and  the  new  vector,  diffmVout  is 
created  for  mod  nt.  Also,  the  vector  named  closeVth/n  shows  the  closest  threshold  level 
for  each  sample  point  of  mod  m.  To  satisfy  the  last  requirement,  the  decimation  region  is 
initially  set  up  for  the  smallest  channel  (mod  9  is  chosen  as  the  smallest  channel  in  our 
project  design).  The  vector  diff9Vout  contains  information  on  the  difference  between  the 
sample  index  and  the  closest  threshold  value.  For  the  upper  parity  threshold  level,  we 
reject  the  sample  if  that  difference  is  smaller  than  the  difference  between  the  original 
threshold  level  and  the  upper  parity  threshold  level.  For  the  lower  parity  threshold  level, 
we  reject  the  sample  if  that  difference  is  smaller  than  the  difference  between  the  original 
threshold  level  and  the  lower  parity  threshold  level.  In  all  other  cases,  we  accept  the 
sample  output.  If  the  sample  is  in  the  rejection  region,  the  parity  bit  shows  a  “1”  (even 
parity),  and  if  the  sample  is  in  the  acceptance  region,  the  parity  bit  shows  a  “0”  (odd 
parity).  Finally,  this  information  is  stored  in  a  file  named  ‘diary’.  Using  the  file, 
sys_simu.m,  we  can  simulate  the  complete  8-bit  SNS  ADC  system,  and  we  can  compare 
the  actual  hardware  performance  using  the  LabVIEW  generated  input  signals.  This  testing 
and  comparison  is  described  in  Chapter  IV. 
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ni.  PARITY  TIMING  PROCESSOR 


A,  OVERVIEW 

In  the  SNS  ADC  hardware,  there  is  a  considerable  difference  in  the  arrival  time  of 
the  parity  bit  and  the  digital  word  representing  the  analog  input.  To  engineer  the  path 
length  to  be  the  same,  a  parity  timing  processor  is  constracted. 

B.  DEVELOPMENT  OF  THE  TIMING  PROCESSOR 

Figure  16  shows  the  system  diagram  of  the  SNS  ADC  that  includes  the  parity 
circuit  information  path,  as  well  as  the  timing  pixx^essor  diagram.  The  timing  processor  is 
divided  into  two  parts:  the  clock  circuit  and  the  holding  circuit.  The  timing  processor  uses 
the  same  detected  laser  pulse  as  an  input  signal.  The  AND  gate  creates  the  clock  signal  for 
the  D-Q  latching  circuit.  The  parity  circuit  signal  is  designed  such  that  an  even  parity  (bit 
information  of  “0”)  gives  an  “ON”,  and  odd  parity  (bit  information  of  “1”)  gives  an  “OFF*. 
So  when  both  the  clock  pulse  and  parity  circuit  signal  are  “ON”,  the  clock  signal  is 
generated  and  drives  the  D-Q  latch.  Otherwise,  the  previous  output  is  held  until  the  next 
clock  pulse.  Figure  17  shows  the  schematic  diagram  of  the  timing  processor.  The  timing 
processor  operates  as  follows:  The  regulator  chip  LM7805CT  generates  the  stable  +5  Vdc 
from  +12  Vdc  to  supply  Vcc  for  all  IC  chips.  The  comparator  AD  9698  cranpares  the  pulse 
input  from  the  photo  detector  with  the  threshold  level  voltage,  and  gives  a  TTL  level  pulse 
output.  The  LM324  is  used  as  a  buffer  for  the  threshold  voltage,  which  is  generated  by  a 
10  kG  potentiometer  and  a  1(X)  G  resistance.  The  PWC-30  accepts  the  TTL  level  pulse  and 
generates  an  output  pulse  of  set  duration.  74S08  is  an  AND  logic  gate  which  generates  the 
clock  pulse  to  drive  the  D-flip  flop  74ALS374.  The  74ALS374  latches  the  digital  word 
from  the  EPROM.  The  74ALS374  is  the  final  stage  of  the  SNS  ADC.  The  following 
subsection  describes  each  device’s  role  in  the  timing  processor  and  specifications. 

L  Laser  Transmitter,  Model  400 

The  laser  transmitter  should  have  as  much  power  as  practical  to  supply  the  system 
with  sufficient  transmit  power.  From  the  previous  research  done  by  LT.  Craig  Crowe,  the 
model  400,  manufactured  by  Broadband  Communication  Products,  Inc.,  is  chosen  as  a 


37 
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ladder 
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Figure  16.  Modified  Block  Diagram  of  the  SNS  ADC 
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Figure  17.  Schematic  Diagram  of  the  Parity  Timing  Processor 
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laser  transmitter  for  this  design.  [Ref.  17]  The  model  400  is  a  high  speed  optical  signal 
source  with  a  data  rate  as  high  as  1.3  giga  bit  per  second  (GBPS).  The  model  400  can  also 
handle  both  digital  and  analog  inputs  up  to  1  GHz.  In  the  digital  mode,  the  model  4(X) 
accepts  an  input  pulse  pattern  from  DC  to  sub-nanosecond  pulses  in  continuous  or  burst 
mode,  and  the  pulse  width  of  the  optical  output  and  las«-  bias  point  can  be  adjusted.  Table 
3  shows  the  characteristic  table  for  the  model  400.  [Ref.  18]  The  model  400  provides  the 
q)ecified  data  rate  and  pulse  width  required  for  this  application.  LT.  Crowe  also  mentioned 
that  the  mode  lock  laser  would  be  a  better  choice  for  the  laser  transmitter  because  it  offers  a 
significant  inq)rovement  in  power  over  that  of  the  model  400.  The  laboratory’s  capability, 
however,  requires  us  to  use  the  model  400.[Ref.  18]  The  optical  pulse  generated  by  the 
model  400  goes  to  a  1  x  4  splitter.  Three  of  the  pulses  go  to  the  folding  circuits,  and  the 
remaining  one  is  detected  and  goes  to  the  timing  processor. 

2.  Photo  Detector,  Model  1811 

The  photo  detector  is  an  optical  to  electrical  converter,  or  optical  receiver.  For  each 
folding  circuit  modulator,  the  wideband  photo  detector  model  310A  was  selected  because  it 
had  the  higher  gain  and  bandwidth  of  any  of  the  equipment  Lt.  Crowe  investigated.  This 
receiver  includes  a  variable  gain.  These  features  are  important  for  the  folding  circuits,  but 
are  not  required  for  the  timing  processor  design.  A  certain  level  of  gain  is  required  in  this 
case  (but  it  does  not  have  to  be  same  gain  as  model  310A):  We  also  must  consider  the 
limitation  of  the  laboratory’s  capability.  For  this  reason,  model  1811,  manufactured  by 
New  Focus,  Inc.,  is  chosen.  The  model  1811  is  a  low-noise  photo  detector  which  uses  an 
InGaAs  PIN  photo  diode  in  the  near  IR  wavelength  region  (800  nm  to  1800  nm).  This  is  a 
range  which  is  includes  the  wavelength  of  the  laser  transmitter.  In  addition,  the  model 
1811  is  more  compact  than  the  model  310A.  The  specifications  of  model  1811  are  shown 
in  Table  4.  [Ref.  19]  The  photo  detector  detects  the  optical  signal  and  converts  it  to  an 
electrical  signal  which  passes  to  the  timing  processor. 

3.  Fast  Speed  Comparator,  AD9698 

When  constructing  an  analog-to-digital  converter,  the  performance  of  the 
comparator  plays  the  key  role  in  the  total  system.  The  AD9698,  manufactured  by  Analog 
Devices,  Inc.,  is  designed  as  an  ultra  fast  TTL-compatible  voltage  comparator.  This  chip 
can  achieve  propagation  delays  previously  possible  only  in  high  performance  emitter- 
coupled  logic  (ECL)  devices.  ECL’s  offer  faster  performance  than  TTL,  but  have 
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Parameter 


Characteristic 


Digital  Input  (Electrical) 


Digital  Output  (Optical) 


Analog  input(Electrical) 

Frequency 

0.1-1,000  MHz 

Input  level 

1.0  V  peak  to  peak  (typ.) 

Input  impedance 

50Q,  AC  coupled 

Analog  Output  (Optical) 

Average  Optical  Power 

0.5  mW 

Frequency  Response 

0.1-1,000  MHz 

Polarity 

Inverting 

Table  3.  Specifications  of  the  BCP  model  400 
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0.5  nsec  (max.) 

0.75  nsec 
10 : 1  (min.) 

0.75  mW  into  (8/125)  fiber 

13(X)nm 

4nm 


Pulse  Rise/Fall  time 
Minimum  Pulse  Width 
Extinction  Ratio 
Peak  Coupled  Power 
Wavelength 
Spectral  Width 


1.3  GBPS 

No  Constraint,  any  combination 


Bit  Rate 
Pulse  Pattern 


1  Parameter 

Characteristic  | 

1  Coupling 

DC  or  AC  1 

Bandwidth  (3  dB) 

DC-125  MHz  (typ.  DC) 

25  kHz-125  MHz  (typ.  AC) 

Wavelength  Range 

800-1800  nm 

Photodiode  Material 

InGaAs  PIN 

Power  Requirements 

±15  Vdc;  250  mA 

Rise  Time 

3  nsec  (typ.) 

Current  Gain 

40V/mA(typ.) 

Input  Noise  Current 

2pA/Hz(l/2)@  10  MHz 

Ouq)ut  Current 

40  mA  (max  into  50  Q) 

Maximum  Input  Power 

5  mW  (max  @1.3  pm) 

Table  4.  Specifications  of  the  model  1811  (New  Focus,  Inc.) 
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problems  with  the  transmission-line  effect  due  to  the  length  and  inductance  of  the  signal 
lines.  Table  5  shows  the  switching  performance  of  the  AD9698.  Rise  time,  ts,  is  one  of 
popular  specifications  used  to  show  the  response  speed  of  the  comparator.  Previous  TTL- 
compatible  fast  comparators  such  as  LT1016  (manufactured  by  Linear  Technology  Corp.) 
or  Am686  (manufactured  by  Advanced  Micro  Devices)  have  6  to  10  nsec  of  tj,  while 
AD9698  has  1.85  nsec  (typ.).  [Ref.  1;  pp.  584-585,  Ref.  20]  Table  6  shows  the 
comparison  of  the  TTL-compatible  comparator  with  other  comparators.  In  the  high  speed 
8-bit  SNS  ADC  design,  many  AD9698  comparators  are  used  in  the  comparator  ladder 
circuit.  In  the  timing  processor,  the  AD9698  is  used  as  the  threshold  detector.  The 
operation  of  the  pulse  detector  is  shown  in  Fig.  17.  The  transmitted  laser  pulse  generated 
by  the  BCP  model  400  comes  to  the  Photo  Detector,  New  Focus  model  1811,  which 
converts  the  optical  signal  to  an  electrical  pulse.  The  amplitude  of  this  pulse  is  small 
compared  to  the  TTL  driving  level.  The  threshold  value  of  the  comparator  is  set  as  a  small 
value  (several  tens  of  mV).  The  comparator  compares  these  small  values  and  creates  an 
ou^ut  signal  at  the  TTL  compatible  level.  For  setting  the  threshold  level,  a  10  turn  -10  kQ 
potentiometer  and  a  100  D  resistor  are  used  to  vary  the  comparator’s  threshold  level  in  the 
mV  region.  The  comparator  is  the  front  end  of  the  timing  processor  with  its  output  going 
to  the  pulse  width  controller  PWC-30. 

4.  Pulse  Width  Controller,  PWC’30 

The  TTL  level  comparator  output  is  sent  to  the  pulse  width  controller,  PWC-30 
series,  manufactured  by  Data  Delay  Device,  Inc.  During  the  process  of  building  up  the 
timing  processor,  it  was  found  that  the  AD9698  has  a  TTL  level  output  representative  of  a 
spike.  This  waveform  causes  the  timing  processor  to  mistrigger.  The  short  pulse  width  of 
the  signal  does  not  have  enough  time  to  hold  the  TTL  level  output  “HIGH”  till  the 
following  stage  starts  operating.  If  the  pulse  cannot  hold  its  ouQ)ut  level  on  HIGH  for  at 
least  the  rise  time  of  the  next  device,  the  system  will  have  no  response.  In  order  to  avoid 
these  problems,  a  pulse  width  controller  is  needed.  Table  7  shows  the  specifications  of  the 
selected  PWC-30  series.  The  features  of  the  PWC-30  are  as  follows: 

•  Exact  control  of  pulse  width  (min.  5  nsec  to  max.  500  nsec) 

•  Rising  edge  trigger 

•  Low  power  consumption 

•  Auto-insertable 
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Parameter 

Testing  temperature 

Characteristics  | 

Propagation  Delay  | 

Input  to  Output  HIGH 

Input  to  Ouq)ut  LOW 

0to70"C 

4.5  nsec  (typ.) 

Latch  Enable  Input  to  Output  HIGH 
Latch  Enable  Input  to  Ou^ut  LOW 

0to+70°C 

6.5  nsec  (typ.) 

Delta  Delay  Between  Ou^uts 

+25  -C 

0.5  nsec  (typ.) 

Propagation  Delay  Dispersion 

20  mV  to  100  mV  Overdrive 

100  mV  to  1.0  V  Overdrive 

+25  *C 

0.1  nsec  (typ.) 

Risellme 

+25  ‘C 

1.85  nsec  (typ.) 

Fall  Time 

+25  *C 

1.35  nsec  (typ.) 

Latch  Enable  | 

Pulse  Width 

+25  X 

2.5  nsec  (typ.) 

Setup  Time 

+25  X 

1.7  nsec  (typ.) 

Hold  Time 

+25  X 

1.9  nsec  (typ.) 

Table  5.  Switching  Performance  of  AD9698  (Analog  Devices,  Inc.) 
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1  Chip  type 

Am686 

LT1016 

EL2019C 

1  Quantity  per 
package 

1 

1 

1 

AD9698  1 

2 

rise-time  (typ.) 

9  (nsec) 

10  (nsec) 

6  (nsec) 

1.85  (nsec) 

Does  the  chip 
cover  Q  and  Q’ 
outyut  ? 

no 

yes 

no 

yes 

Is  the  chip 

O.K.  for 
supplying 
single  5  V  ? 

no 

yes 

no 

yes 

Table  6.  Comparison  of  the  Fast  TTL  Level  Comparator 
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Comparator 


Figure  18.  Operation  of  the  Threshold  Detector 
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Parameters 

Characteristics 

Trigger  inherent  delay 

Too  =  6.5  ±  1  nsec 

Pulse- width  tolerance 

5%  or  2  nsec 

Maximum  Input  PRR 

PW  +  20  nsec 

Output  rise-time 

2  nsec  (Typ.) 

HIGH-level  input  voltage 

2.0  V  (min.) 

LOW-level  input  voltage 

+  0.8  V  (max.) 

fflGH-level  output  voltage 

3.4  V  (typ.),  2.7  V  (min.) 

LOW-level  output  voltage 

+0.5  V  (max) 

Table  7.  Specifications  of  the  PWC-30  (Data  Delay  Device,  Inc.) 
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•  High  speed 

Because  it  uses  a  rising  edge  trigger,  the  PWC-30  can  avoid  the  multiple  trigger  problem. 
Figure  19  shows  the  block  diagram  of  the  PWC-30.  In  this  timing  processor,  the  PWC- 
30-50,  generates  50  nsec  of  pulse  width. 

5.  Programmable  Delay  Units,  PDU~16F 

This  device  is  an  optional  device  in  the  timing  processor.  At  the  final  stage  of  the 
SNS  ADC,  the  D-Q  latch  should  handle  both  the  digital  word  information  from  the 
EPROM  and  the  clock  signal.  The  data  path  which  handles  the  digital  word  and  the  path  of 
the  timing  processor  signal  do  not  have  same  propagation  delay.  The  EPROM  has  the 
largest  propagation  delay  (around  1(X)  nsec)  in  the  data  path,  while  the  timing  processor  has 
much  less  of  a  propagation  delay.  In  order  to  get  both  data  and  timing  information  at  the 
same  time,  we  need  to  insert  delay  devices.  There  are  many  options  to  set  up  a  delay  lime. 
An  easy  way  is  to  use  an  optical  delay  line  (ODL).  In  this  method,  an  additional  length  of 
fiber  cable  is  inserted  in  the  optical  signal  line  to  get  the  desired  time  delay.  The  ODL, 
howevCT,  offers  little  flexibility  in  tuning  the  delay  time.  One  ODL  corresponds  only  to  one 
fixed  delay  time.  During  the  process  of  building  up  the  timing  processor,  we  cannot  expect 
to  know  the  exact  time  delay  of  both  the  data  and  timing  signal.  Using  an  electrical  delay 
unit  is  another  way  of  creating  a  specific  delay  time.  Data  Delay  Devices,  Inc.,  offers 
various  types  of  delay  units.  Programmable  delay  units  (PDU)  were  chosen  as  the 
electrical  delay  line  device.  The  PDU  series  is  a  fast  logic  programmable  delay  unit.  The 
PDU  has  a  3-bit  to  8-bit  programmable  delay  line,  and  we  can  choose  the  suitable  chip 
among  the  various  incremental  delay  step  sizes.  In  this  timing  processor,  PDU-16F,  which 
has  6-bit  programmable  delay  lines,  is  used.  The  block  diagram  of  PDU-16F  is  shown  in 
Fig.  20.  For  the  timing  processor,  we  prepare  PDU-16F-2,  which  has  2  +  0.5  nsec  of 
incremental  delay  for  each  step  and  a  total  programmed  delay  of  126  nsec.  The 
specifications  of  the  PDU  are  shown  in  Table  8. 
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Figure  19.  Block  Diagram  of  the  PWC-30  (Data  Delay  Devices,  Inc.) 
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Figure  20.  Block  Diagram  of  the  PDU-16F  (Data  Delay  Devices,  Inc.) 
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Parameter 

Characteristics 

Propagation  Delay 

Address  to  ou^ut  (Tsua) 

7  nsec  (typ.) 

Enable  to  output  (Tsue) 

6  nsec  (typ.) 

Delay  variation 

Monotonic  in  one  direction 

Programmed  delay  tolerance 

±  5%  or  1  nsec 

Inherent  delay 

9  nsec  (pin  number  2)  (typ.) 

8  nsec  (pin  number  1)  (typ.) 

Minimum  pulse- width 

10%  of  total  delay 

HIGH-level  ouq)ut  voltage 

2.5  V(min.)  3.4  V  (typ.) 

LOW-level  ou^ut  voltage 

0.5  V(max.)  0.35  V(typ.) 

HIGH-level  input  voltage 

2.0  V  (min.) 

LOW-level  input  voltage 

0.8  V  (max.) 

Table  8.  Specifications  of  the  PDU-16F  (Data  Delay  Devices,  Inc.) 
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IV.  RESULTS 


A.  8-BIT  SNS  ADC  PERFORMANCE 

L  Overview  of  Decimation 

The  original  SNS  ADC  can  present  a  problem  at  each  of  the  code  transition  points. 
The  probability  of  this  occurrence  depends  upon  the  offset  voltage  match  of  the 
comparators,  the  tolerance  of  the  voltage  dividing  resistors,  and  the  corresponding  width  of 
the  input  voltage  region  being  affected.  The  modified  SNS  ADC  with  a  parity  circuit,  is 
introduced  to  solve  this  problem.  The  parity  circuit  is  used  to  determine  if  the  sample  lies  at 
a  code  transition  point.  The  goal  of  this  section  is  to  test  the  modified  SNS  ADC  hardware 
using  the  LabVEBW-genCTated  input.  Also,  a  comparison  of  the  actual  circuit  response  to 
the  simulation  is  investigated.  Several  decimation  widths  are  tested  to  determine  the 
effectiveness  of  the  parity  circuit. 

2.  Testing  Configuration 

The  system  configuration  for  testing  the  circuit  is  shown  in  Fig.  21.  We  use  an  8- 
bit  SNS  ADC  with  the  LM311  general  purpose  comparator,  which  has  a  200  nsec  (typ.) 
response  time.  The  SNS  ADC  uses  a  “static”  input  sampling  signal  generated  by  the 
LabVIEW  program,  CRAIG7.VI.  For  a  “static”  signal,  each  input  voltage  level  has  a 
certain  time  period  which  is  relatively  long.  For  CRAIG7.VI,  we  choose  a  time  step  of  0.5 
sec  with  the  LabVIEW  signal  output  changing  levels  at  each  new  time  step.  The  system 
moduli  for  the  SNS  ADC  are  mod  9,  10,  and  11.  The  corresponding  total  number  of 
conqjarators  is  37  with  a  maximum  of  18  loaded  in  parallel.  The  mod  9  channel  is  used  for 
the  parity  circuit,  and  the  chosen  decimation  widths  of  0%  (no  decimation),  10%,  15%, 
20%,  and_30%  were  investigated  using  the  corresponding  resistance  networks.  Values  for 
the  folding  circuit  resistance  ladder  (except  for  the  parity  circuit),  and  the  target  threshold 
voltages,  are  shown  in  Table  9.  The  values  for  the  parity  circuit  resistance  ladder  and  target 
threshold  voltages  are  shown  in  Table  10.  No  timing  circuit  was  required  for  these  tests, 
since  the  LabVIEW  program  controls  the  input  voltage  sequence,  and  this  is  a  static 
analysis  of  the  system.  The  signal  flow  of  the  testing  is  as  follows:  First,  CRAIG7.VI 
creates  the  input  folding  waveforms  from  the  SNS  ADC  with  a  0.5  sec  time  step  between 


53 


Figure  21.  Block  Diagram  of  Testing  8-bit  SNS  ADC  Circuit 


mod  10  mod  11 


Reference  Voltage 

5.00  (V) 

5.00  (V) 

Signal  Input  Range 

0to5.0(V) 

0to5.0(V) 

Ladder  resistance  (kf2) 

total  resistance  Rtot 

=  100  kfl  for  modulus 

Resistance  No.l 

2.45 

2.03 

Resistance  No.2 

7.10 

5.91 

Resistance  No.3 

11.10 

9.32 

Resistance  No.4 

11.90 

11.97 

Resistance  No.5 

15.50 

13.66 

Resistance  No.  6 

15.50 

14.23 

Resistance  No.7 

11.90 

13.66 

Resistance  No.  8 

11.10 

11.97 

Resistance  No.9 

7.10 

9.32 

Resistance  No.  10 

2.45 

5.91 

ResistanceNo.il 

— 

2.03 

Table  9.  Testing  Configuration  of  Folding  Circuit  (Value  of  Resistance  Ladder) 
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mod  10 

mod  11  1 

1  Target  Threshold  (V)  | 

Threshold  Voltage  No.l 

0.123 

0.102 

Threshold  Voltage  No.2 

0.478 

0.397 

Threshold  Voltage  No.3 

1.031 

0.863 

Threshold  Voltage  No.4 

1.728 

1.462 

Threshold  Voltage  No.5 

2.500 

2.144 

Threshold  Voltage  No.6 

3.273 

2.856 

Threshold  Voltage  No.7 

3.970 

3.538 

Threshold  Voltage  No.8 

4.523 

4.137 

1  Threshold  Voltage  No.9 

4.878 

4.603 

1  Threshold  Voltage  No.  10 

-- 

4.899 

Table  9.(cont.)  Testing  Configuration  of  Folding  Circuit  (Threshold  Voltage ) 
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mod9 

(Parity 

circuit) 

0% 

decimation 

10% 

decimation 

15% 

decimation 

20% 

decimation 

30% 

decimation 

1  Resistor  Rtot  =  100  kQ  for  each  case  | 

R1 

— 

10 

17 

30 

R2 

3,020 

2,716 

2,567 

2,417 

2,116 

R3 

— 

597 

895 

1,194 

1,790 

R4 

— 

8,680 

_ 

7,822 

7,390 

6,958 

6,094 

R5 

— 

1,122 

1,683 

2,244 

3,364 

R6 

13,300 

11,983 

11,323 

10,661 

9,336 

R7 

— 

1,511 

2,267 

3,022 

4,532 

R8 

16,320 

14,700 

13,890 

13,078 

11,452 

R9 

— 

1,719 

2,578 

3,437 

5,154 

RIO 

17,360 

15,643 

13,890 

13,078 

11,452 

Rll 

— 

1,719 

2,578 

3,437 

5,154 

R12 

16,320 

14,700 

13,890 

13,078 

11,452 

R13 

— 

1,511 

2,267 

3,022 

4,532 

R14 

13,300 

11,983 

11,323 

10,661 

9,336 

R15 

— 

1,122 

1,683 

2,244 

3,364 

R16 

8,680 

7,822 

7,390 

6958 

6.094| 

R17 

— 

597 

895 

1,194 

l,79o| 

R18- 

3,020 

2,716 

2,567 

2,417 

2,1161 

R19 

10 

17 

30 

681 

Table  10.  Testing  Configuration  of  Parity  Circuit  (Value  of  Resistance  Ladder) 
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mod9 

(Parity  decimation 
circuit) 


Target  Threshold 


Vthl 


0.151  (V) 


Signal 

Input 

Range 


0.38  (mV) 

0.86  (mV) 

136.20 

129.19 

166.05 

173.96 

557.14 

543.48 

613.23 

627.61 

1,212.40 

1,193.75 

1,287.98 

1,307.10 

2.023  (V) 

2.002  (V) 

2.109 

2.130 

2.891 

2.870 

2.977 

2.998 

3.712 

3.693 

3.786 

4.387 

4.372 

4.443 

4.457 

4.834 

4.826 

4.864 

4.871 

4.999 

4.999 

5.00(V) 

5.00(V) 

0  to  5.0(V) 

0to5.0(V) 

i 

20% 

decimation 

30% 

decimation 

1.52  (mV) 

3.43  (mV) 

122.36 

109.24 

198.74 

503.41 

642.14 

671.66 

1,175.20 

1,138.40 

1,326.32 

1,365.02 

1,980.22 

1,937.62 

2.152  (V) 

2.195  (V) 

2.848 

2.805 

3.018 

3.062 

3.674 

3.635 

3.825 

3.862 

4.358 

4.328 

4.470 

4.497 

4.818 

o 

oo 

4.878 

4.891 

4.998 

4.997 

5.00(V) 

5.00(V) 

0to5.0(V) 

0to5.0(V) 

Table  10.(cont.)  Testing  Configulation  of  Parity  Circuit  (Threshold  Voltage) 


samples.  The  SNS  ADC  converts  that  signal  to  a  thermometer  code  using  the  comparator 
ladder  and  changes  the  thermometer  code  to  a  binary  code  by  means  of  a  program  logic 
array  (PLA)  such  as  an  EPROM.  The  mod  9  comparator  lines  are  also  used  to  generate  the 
parity.  To  process  this  information,  a  Hewlett  Packard  model  1631  AID  logic  analyzer  is 
used.  The  model  1631  A/D  is  a  general  purpose  logic  analyzer  for  performing  state,  timing, 
and  analog  waveform  measurements.  The  output  data  and  parity  bit  information  are  stored 
in  the  storage  device  (floppy  disk)  via  the  logic  analyzer.  The  data  for  each  decimation 
width  is  stored  and  compared.  This  information  is  also  compared  to  the  output  of  the 
MATLAB  simulation  program,  sys_simu.m. 

3.  Testing  Results 

CRAIG7.VI  creates  a  cosine  square  wave  with  a  5.0  V  peak  to  peak  voltage 
amplitude.  Each  target  threshold  voltage  and  resistance  value  is  computed  by  MATLAB 
(Vth_calc.m,  and  idlprty.m).  A  combination  of  fixed  resistors  and  potentiometers  are  used 
in  the  con5)arator  resistance  ladder  to  adjust  the  threshold  level.  The  total  resistance  value 
of  each  ladder  is  100  kD.  The  hardest  part  of  building  up  the  circuit  is  adjusting  these 
threshold  voltages.  Adjusting  one  element  in  the  ladder  circuit  to  get  a  desired  threshold 
voltage  affects  all  other  threshold  values,  so  the  remaining  resistances  need  to  be  adjusted 
as  well. 

Adjusting  threshold  voltages  for  the  parity  circuit  modulus  were  especially  difficult 
since  the  circuit  has  pairs  of  resistance  values  with  only  a  small  difference.  In  addition,  it 
has  2m  1  potentiometers  for  the  smallest  channel  of  mod  mi.  The  logic  analyzer  is  used  to 
check  all  of  the  logic  values  from  the  SNS  ADC.  Figure  22  shows  the  comparator  ladder 
output  of  mod  10  captured  from  the  logic  analyzer  display.  Clearly  the  output  is 
symmetric.  Labels  BO  to  B8  of  the  figure  indicate  each  comparator  output.  BO  shows  the 
bottom  (or  lowest)  level  comparator  state,  and  B8  shows  the  top  (or  highest)  level 
comparator  state.  After  adjusting  the  resistance  value  to  get  the  symmetrical  output,  the 
static  analysis  is  performed.  The  output  state,  decimal  number  representation  and  parity  bit 
information,  are  collected  for  each  decimation  width  (10%,  15%,  20%,  and  30%).  Table 
1 1  shows  the  parity  circuit  behavior  of  the  8-bit  SNS  ADC.  MATLAB  is  used  as  a 
reference  for  the  parity  bit  information.  The  parity  “1”  bit  indicates  a  rejected  sample 
proportional  to  the  decimation  width.  In  this  test,  646  samples  are  used  to  express  the  8- 
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Table  11.  Parity  Circuit,  Decimation  and  Hardware  Performance  Comparison 


bits.  As  the  decimation  width  increases,  the  total  number  of  parity  “1”  bits  increases,  and 
the  recognition  rate,  which  shows  how  much  the  actual  number  of  parity  “1”  bits  agree 
with  the  MATLAB-predicted  parity  “1”  bits,  also  increases.  The  larger  the  decimation 
width,  the  fewer  number  of  output  errors.  However,  the  larger  the  decimation  width,  the 
higher  the  number  of  samples  rejected.  For  example,  in  Table  11, 30%  decimation  gives  a 
recognition  rate  in  excess  of  95%,  but  over  one  third  of  the  total  samples  are  rejected. 
From  Table  11,  it  can  be  said  the  rate  of  rejected  samples  versus  the  total  number  of 
samples  remains  almost  constant  for  decimation  widths  less  than  20%,  and  recognition 
rates  remain  constant  till  the  15%  decimation  width. 

The  MATLAB  file  errormod.m  is  used  as  a  comparison  tool  for  “before”  and 
“after”  corrections.  The  logic  analyzer  model  1631  A/D  stores  the  sample  index,  output 
state  (decimal),  and  the  parity  bit  data.  The  file  errormod.m  checks  these  three  sets  of 
information,  and  it  regards  the  data  as  corrected  when  the  output  state  error  and  parity  bit 
“1”  occur  at  the  same  time.  If  the  output  error  occurred,  and  the  parity  bit  was  “0”,  then 
MATLAB  regards  the  error  as  uncorrected.  This  means  the  system  put  out  erroneous 
information.  Obviously  a  higher  error  correction  rate  is  better.  Table  12  shows  the 
e3q)erimental  results  of  error  correction  rates.  Each  decimation  width  case  is  run  three  times 
to  collect  data.  The  larger  decimation  width  correlates  with  a  better  eiror  correction  rate. 

Figures  23  through  26  show  the  results  of  error  correction  for  each  decimation 
width.  For  each  case,  the  errors  present  decrease  significantly,  and  a  20%  decimation 
width  is  enough  to  get  a  100%  error  correction.  In  these  figures,  the  x-axis  shows  the 
sample  index  which  relates  to  the  input  voltage  of  each  MZI,  and  the  y-axis  shows  the 
decimal  output  state.  Some  of  the  errors  are  within  the  8-bit  state  (less  than  256).  Others 
are  out  of  the  8-bit  state  (greater  than  or  equal  to  256).  Figure  27  shows  the  relationship  of 
error  correction  and  decimation  width  graphically.  The  results  are  summarized  as  follows: 

•  -  The  larger  the  decimation  width  the  higher  the  recognition  rate  of  the 

expected  parity  bit 

•  Larger  decimation  widths  yield  higher  error  correction  rates. 

•  Larger  decimation  widths  may  cause  excessive  problems  such  as  high  rates 
of  rejecting  sample  output. 

20%  decimation  width  is  recommended  as  the  best  choice  for  the  parity  circuit 
configuration. 
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%Decimation 

Total  number  of 

OTor 

Total  numbCT  of 

corrected  error 

Error  ccMTection 

rates  (%) 

Average  error 

correction  rates 

10 

18 

16 

88.8 

84.1 

10 

21 

17 

80.1 

10 

24 

20 

83.3 

15 

18 

18 

100 

89.8 

15 

24 

23 

95.8 

15 

18 

16 

88.8 

20 

28 

18 

100 

100 

20 

37 

18 

100 

20 

30 

22 

100 

30 

23 

23 

100 

100 

30 

19 

19 

100 

30 

23 

23 

100 

Table  12.  Error  Correction  Data  of  8-bit  SNS  ADC 
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output  decimal  number  output  decimal  number 


Figure  23.  Characteristic  Analysis  of  an  8-bit  SNS  ADC  (10%Decimation) 
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output  decimal  number  output 


Figure  24.  Characteristic  Analysis  of  an  8-bit  SNS  ADC  (15%Decimation) 
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output  decimal  number  output  decima!  number 


Figure  25.  Characteristic  Analysis  of  an  8-bit  SNS  ADC  (20%Decimation) 
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output  decimal  number  output  decimal  number 


Input  voltage  index 


Figure  26.  Characteristic  Analysis  of  an  8-bit  SNS  ADC  (30%Decimation) 
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Error  U 


B.  TESTING  AND  EVALUATING  THE  TIMING  PROCESSOR 


L  Signal  Performance 

The  timing  prcx:essor  is  built  as  discussed  in  Chapter  m.  A  picture  of  the  timing 
circuit  is  shown  in  Fig.  28.  The  circuit  is  made  by  using  a  wire  wrapping  board.  Supply 
voltages  (+12  V,  ±5  V)  and  ground  for  the  circuit  are  applied.  The  input  signal  and  the 
output  of  each  stage  can  be  detected  from  the  BNC  terminal  in  this  circuit.  Two  DIP 
switohes  are  shown  on  the  circuit.  One  is  for  simulating  the  EPROM  output  signal,  and  the 
other  is  used  as  an  address  line  for  the  programmable  delay  units,  PDU-16F.  In  order  to 
cool  down  the  chip,  the  regulator  LM7805CT  has  a  heat  sink.  Capacitors  are  added  for 
each  chip  to  terminate  the  high  frequency  effects.  A  Digitizing  Oscilloscope  TDS420  is 
used  to  capture  the  waveform  of  the  timing  processor.  The  TDS420  can  take  snapshots  of 
the  four  individual  waveforms.  Each  terminal  can  choose  a  50  Q  or  1  MO  terminal.  The 
TDS  420  stores  the  image  in  the  personal  computer  via  a  GPIB  bus.  Also  the  TDS420  can 
handle  a  100  mega  sample  per  second  (MSPS)  maximum  digitizing  rate,  and  a  150  MHz 
maximum  analog  bandwidth.  [Ref.  22]  Figure  29  shows  the  output  waveform  of  the  photo 
detector,  model  1811.  The  laser  transmitter  model  400  applies  a  20  nsec  pulse  with  a  pulse 
repetition  interval  (PRI)  of  200  nsec  in  this  figure  (10%  duty  cycle).  The  figure  shows  that 
the  detected  pulse  width  is  almost  the  same  as  the  pulse  width  of  the  laser  transmitter, 
which  means  the  propagation  delay  is  minimal.  The  figure  also  shows  that  the  highest 
value  of  the  photo  detector  is  146.0  mV,  and  the  average  noise  layer  value  is  around  15 
mV.  This  input  signal,  however,  changes  its  waveform  when  the  photo  detector  output 
connects  to  the  input  terminal  of  the  timing  processor.  Figure  30  shows  the  signal 
waveform  at  the  input  terminal  of  the  timing  processor  without  applying  power.  Instead  of 
the  nice  pulse  shape  of  Fig.  29,  there  are  some  higher  order  frequency  effects  in  the 
waveform.  This  waveform  also  changes  when  supply  voltage  is  applied.  Figure  31  shows 
the  signal  waveform  at  the  input  tominal  of  the  timing  processor  after  supplying  the  power. 
The  waveform  in  Fig.  31  clearly  shows  multiple  peaks.  The  maximum  measured  voltage 
of  the  main  pulse  is  360  mV,  but  some  of  the  other  pulses  have  voltage  spikes  in  excess  of 
200  mV.  This  input  signal  passes  through  the  AD9698  comparator.  However,  it  is  hard  to 
set  up  the  threshold  voltage  of  the  comparator  to  get  the  main  pulse  only,  since  the  main 
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Figure  28.  Picture  of  the  Timing  Processor 
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Figure  29.  Output  Waveform  of  the  Photo  Receiver  model  1811(New  Focus,  Inc.) 
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Figure  30.  Input  Signal  Waveform  of  the  Timing  Processor  Without  Supply 
Voltage 
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pulse  has  a  dip  whose  value  is  close  to  the  peak  value  of  the  second  pulse.  This  causes 
multiple  triggering  of  the  conparator,  and  consequently  a  plural  TTL  level  ouq)ut  from  the 
comparator.  The  ouq)ut  waveform  of  the  comparator  AD9698  is  shown  in  Fig.  32.  The 
AD9698  generates  three  sharp  pulses  because  of  the  multiple  triggers  as  shown  in  Fig.  31. 
These  pulses  have  a  TTL  level  output  and  also  are  very  sharp.  Each  pulse  has  a  5  to  6  nsec 
pulse  width,  so  this  usually  causes  a  misfire  of  the  clock  and  a  big  timing  problem.  This 
effect,  however,  is  diminished  by  making  use  of  the  pulse  width  controller  PWC-30.  As 
noted  in  ChaptCT  IQ,  the  PWC-30  has  a  very  effective  feature  in  that  it  activates  using  a 
rising  edge  trigger,  and  no  other  pulse  can  retrigger  the  device  during  the  operation.  Figure 
33  shows  the  charactaistics  of  the  PWC-30  in  operation.  This  figure  explains  how  the 
PWC-30  works  for  various  inputs.  The  signal  output  of  the  PWC-30  is  independent  of  the 
pulse  width  of  the  input  signal  and  each  signal  output  has  a  fixed  pulse  width  x.  The 
waveform  difference  of  the  PWC-30  input  /  output  is  shown  in  Fig.  34.  The  upper  signal 
is  the  input  to  the  PWC-30  and  this  waveform  is  almost  the  same  as  the  one  in  Fig.  32. 
The  lower  waveform  is  the  PWC  output  and  this  pulse  has  a  50  nsec  pulse  width.  By 
measurement,  the  propagation  delay  of  the  PWC-30  is  around  7  nsec  which  is  within  the 
specifications.  The  output  signal  shows  two  dips  in  the  pulse,but  these  local  minima  are 
both  above  the  “ON”  level,  and  therefore  do  not  affect  other  parts  of  the  timing  processor. 
The  50  nsec  pulse  signal  generated  by  the  PWC-30  and  the  parity  circuit  output  both  go 
through  the  AND  logic  (74S08).  Finally,  the  clock  pulse  of  the  D-Q  latch  is  obtained.  The 
timing  processor  and  folding  circuit  could  not  be  integrated  at  this  time,  so  a  simulated 
parity  signal  was  used  in  testing  instead  of  the  actual  parity  signal.  In  the  following 
section,  the  simulated  parity  signal  is  generated  by  WAVETEK  model  145  or  by  toggling 
the  ouq)ut  of  the  PWC-30.  WAVETEK  model  145  can  generate  a  TTL  level  pulse  signal  at 
frequencies  up  to  20  MHz. 

2.~  Timing  Analysis  of  the  Processor 

The  next  step  in  the  evaluation  of  the  timing  processor  is  to  investigate  the  riming 
relationship  for  each  stage  of  the  processor,  and  to  confirm  the  operation  of  the  timing 
processor.  Generally,  each  IC  chip  has  a  certain  propagation  delay  time  during  the  signal 
VO.  Even  the  fastest  devices  experience  these  delays.  In  the  timing  processor,  several 
chips  are  used,  and  it  is  important  to  know  the  total  system  time  delay.  If  the  time  delay  is 
longer  than  the  PRI,  then  the  sampling  is  not  working  correctly.  The  system  does  not 
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Figure  33.  Operation  Characteristics  of  the  PWC-30 
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work  because  the  timing  clock  does  not  match  the  sample.  A  similar  problem  occurs  when 
the  time  delay  is  too  short  compared  to  the  data.  However,  we  want  to  avoid  the  former 
case  rather  than  the  latter  case  because  the  data  bus  has  many  lines.  Qjrrecting  for  the  time 
delay  in  the  latter  case  would  require  only  one  delay  line  for  the  time  delay.  Figure  35 
shows  the  clock  pulse  waveform  of  each  level.  Channel  1  shows  the  input  signal  to  the 
system.  The  rising  edge  of  this  signal  is  used  as  the  reference  point  of  the  propagation  time 
delay  analysis.  Channel  2  shows  the  output  of  the  PWC-30.  Channel  3  has  the  same 
signal  path,  but  it  is  measured  at  the  AND  logic  input.  The  last  channel  shows  the  output 
of  the  AND  logic  which  is  the  actual  clock  signal  of  the  D-Q  latches.  In  this  case,  the 
output  signal  of  the  PWC-30  is  applied  to  both  74S08  input  pins  as  a  timing  clock  pulse 
and  a  simulated  parity  check  signal.  The  right  side  of  the  figure  shows  the  delay  time  of 
each  channel  compared  to  channel  1.  According  to  the  figure,  the  total  delay  time  from 
photo  detector  to  the  final  clock  output  is  50  nsec.  If  we  obtain  the  propagation  time  delay 
of  each  component  firom  the  specification  data,  the  total  delay  time  should  be  around  20 
nsec.  The  timing  delay,  however,  can  be  affected  by  the  wiring.  The  biggest  cause  of  the 
delay  is  the  signal  line  from  the  photo  detector  to  the  tiirdng  processor  since  a  very  long 
coaxial  cable  is  used  in  this  experiment.  Even  the  short  length  of  wire  used  in  the  wire 
wrapping  can  contribute  to  the  total  delay.  For  example,  in  an  ideal  circuit,  the  pulses  in 
chaimels  2  and  3  should  be  coincident  and  have  the  same  timing  delay  compared  to  the 
photo  detector  output.  However,  Fig.  35  shows  an  8.5  nsec  pulse  delay  between  the  two 
devices  in  the  actual  circuit.  The  final  clock  pulse  has  a  TTL  level  output,  and  its  pulse 
width  is  50  nsec.  This  signal  arrives  at  the  clock  input  pin  of  D-flip  flop  74ALS374,  and 
proceeds  to  update  the  data  output.  In  this  configuration,  the  PRI  of  the  system  is  set  at 
200  nsec,  so  the  timing  process  will  eliminate  the  mismatching  of  sampling  data  and  timing 
signal  at  the  D-flip  flop.  This  timing  process  is  also  analyzed  using  the  simulated  parity  bit 
signal,  and  simulated  data  bus  information.  The  Hewlett  Packard  logic  analysis  system 
model  16500B  is  used  for  timing  analysis.  The  model  16500B  has  much  more  capability 
than  the  model  1631A/D  logic  analyzer  used  in  the  decimation  error  testing.  Model  16500B 
has  some  features  not  included  in  the  model  1631  A/D  such  as: 

*  100  MHz  state  and  500  MHz  timing  acquisition  speed 

•  96  data  channels  /  6  clocks 

♦  GPIB  capability  (programming,  data  stored) 
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The  waveform  for  each  channel  indicates  the  following  signal 

Channel  1  :  the  input  signal  to  the  system  (from  the  photo  detector) 
Channel  2  :  the  output  waveform  of  the  PWC-30 
Channel  3  :  the  input  waveform  of  the  74S08  (same  path  of  the  Ch.  2) 
Channel  4  :  the  output  wave  form  of  the  74S08 _ 


Figure  35.  Waveform  of  Each  Stage  of  the  Timing  Clock  Signal 
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♦  3.5  nsec  window 

This  equipment  is  mainly  used  for  SNS  ADC  folding  circuits  analysis.  [Ref.  23]  The 
WAVETEK  model  145  is  used  as  a  simulated  parity  bit  signal.  Figure  36  shows  the 
display  of  the  model  16500B.  The  index  numbers  0  to  3  are  used.  Each  signal  is  a  TTL 
level  output,  and  each  index  shows  the  ouq)ut  of  the  comparator  AD9698,  the  output  of  the 
pulse  width  controller  PWC-30,  the  simulated  parity  bit  signal  from  WAVETEK  model 
145,  and  the  final  clock  pulse  to  the  D-flip  flop.  The  laser  pulse  configuration  is  the  same 
as  the  previous  experiment  with  a  20  nsec  pulse  width  and  a  PRI  is  equal  to  200  nsec.  The 
model  145  generates  a  700  kHz  square  wave  with  a  50%  duty  cycle  to  confirm  the 
operation  of  the  AND  logic.  The  multiple  triggering  of  the  comparator  is  seen  on  the  index 
0  as  mentioned  in  the  previous  section.  The  pulse  width  controller,  however,  eliminates 
this  problem,  and  generates  a  periodic  pulse  train  with  a  50  nsec  pulse  width.  The  74S08 
AND  logic  compares  the  clock  input  from  the  PWC-30  with  parity  bit  information,  and  the 
final  output  is  the  exact  pulse  waveform  that  we  predict.  Another  example  of  the  timing 
analysis  is  shown  in  Fig.  37.  The  model  145  generates  a  5  MHz  pulse  with  a  50%  duty 
cycle  as  a  parity  bit  signal.  Index  5  shows  the  ouq)ut  waveform  of  the  programmable  delay 
imit,  PDU-16F.  The  pulse  delay  unit  maintains  a  certain  time  delay  and  pulse  width  from 
the  74S08  output.  For  the  programmable  delay  unit,  a  DIP  switch  is  used  to  enable  each 
address  line  to  change  the  delay  time.  The  fixed  delay  time  is  explicit  in  the  figiue.  In  this 
example,  address  line  number  5  is  in  the  ON  state,  and  the  OFF  state  signal  is  applied  to 
other  address  lines.  According  to  the  specifications  for  the  programmable  delay  unit  in 
Table  8,  the  inherent  delay  on  pin  1  is  8  nsec,  and  the  incremental  delay  ATd  is  2  ±  0.5  nsec 

per  step.  Table  13  shows  the  tmth  table  of  the  PDU-16F.  According  to  Table  13,  the 
configuration  of  this  example  produces  a  delay  output  Tie.  The  relationship  between  the 
step  number  na  and  the  time  delay  Ta  is 

Tj  —  T D0+  ATj'n^  (26) 

Therefore  the  Ta  of  a  delay  output  Tie  is  40  ±  8  nsec.  By  observation  of  Fig.  37  using  the 
timing  cursor,  the  time  difference  between  two  signals  (input  /  output  of  PDU-16F)  is  44.0 
nsec.  This  satisfies  the  specifications,  and  confirms  that  the  PDU-16F  is  operating 
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Each  label  shows  the  following  TTL  level  output : 

Label  0  :  the  output  signal  of  the  AD9698 

Label  1  :  the  output  signal  of  the  PWC-30 

Label  2  :  the  simulated  parity  bit  signal  (WAVETEK  145) 

Label  3  :  the  output  signal  of  the  74S08 

(Label  4  to  7  :  unconnected) _ 


Figure  36.  Display  of  the  Logic  Analyzer  Model  165006 
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Each  label  shows  the  following  TTL  level  output : 

Label  0  :  the  output  signal  of  the  AD9698 

Label  1  :  the  output  signal  of  the  PWC-30 

Label  2  :  the  simulated  parity  bit  signal  (WAVETEK  145) 

Label  3_:  the  output  signal  of  the  74S08 

Label  4  :  the  output  signal  of  PDU-16F 

(Label  5  to  7  :  unconnected) _ 


Figure  37.  Another  Example  of  the  Timing  Signal  Analysis 
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•  0 :  Ul'F  state  1 :  OJN  state 


*  Don’t  care  condition  if  enable  (Eo)  is  ON 
Table  13.  Truth  Table  of  the  PDU-16F  Address  Line 
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correctly.  Figure  38  shows  the  timing  signal  behavior.  The  propagation  delay  of  each 
stage  in  TTL  level  logic  is  measured  using  the  same  manner  as  in  the  previous  example. 
Figure  37  differs  from  Fig.  36  in  that  the  rising  edge  of  the  parity  bit  information  and  the 
PWC-30  output  are  not  totally  overlapped.  Therefore,  the  74S08  output  has  a  shorter  pulse 
width.  This  is,  however,  closer  to  the  actual  signal  behavior  because  the  actual  parity  bit 
signal  has  a  narrower  pulse  width  than  the  simulated  parity  bit  signal.  Using  the  two 
cinsors  shown  in  the  figure,  we  measured  the  propagation  delay  with  the  results  listed  in 
Table  14.  The  first  row  in  this  table  shows  the  propagation  delay  of  the  PWC-30.  The 
second  row  indicates  the  propagation  delay  of  74S08. 

Finally  the  simulated  EPROM  output  data  is  applied  to  the  74ALS374  octal  D-flip 
flop  and  the  input  /  ouq)ut  of  74ALS374  is  observed.  The  data  output  states  are  controlled 
by  a  dip  switch,  and  this  data  signal  is  static.  Each  state  signal  is  a  DC  output  and  not  a 
high  frequency  signal  ouq)ut.  Figure  38  shows  one  set  of  results.  Index  labels  0  to  7 
show  the  data  input  of  the  74ALS374,  and  index  labels  8  to  15  show  the  data  output  of  the 
74ALS374.  The  DIP  switch  is  set  to  generate  the  binary  word  “10101010”.  The  word  is 
applied  to  the  input  of  the  74ALS374.  The  output  should  be  the  same  binary  word 
“10101010”,  and  the  same  information  does  in  fact  come  out  for  each  corresponding  pin. 
Because  of  the  static  input,  the  data  for  transient  I/O  behavior  in  the  circuit  cannot  be 
recorded,  and  characteristic  analysis  using  optical  pulse  signal  is  needed  for  further 
research. 
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Each  label  shows  the  following  TTL  level  output : 

Label  0  :  the  output  signal  of  the  AD9698 
Label  1  :  the  output  signal  of  the  PWC-30 
Label  2  :  the  simulated  parity  bit  signal  (WAVETEK  145) 

Label  3  :  the  output  signal  of  the  74S08 
Label  4  :  the  output  signal  of  PDU-16F 

(Label  5  to  7  :  unconnected) _ 


Figure  38.  Details  of  the  Timing  Processor  Signal  Behavior 
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Rise  Edge  Reference  A 

Rise  Edge  Reference  B 

Propagation  Delay 

AD9698  output 

PWC-30  output 

8.0  (nsec) 

PWC-30  output 

74S08  output 

4.0  (nsec) 

74S08  output 

PDU-16F  output 

44.0  (nsec)* 

*  Choose  Delay  Out  as  Tje  and  this  satisfies  the  specification 


Table  14.  Observation  Results  of  Propagation  Delay 
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V.  SUMMARY  AND  FURTHER  RESEARCH 


This  thesis  provides  details  on  the  system  simulation  of  an  8-bit  SNS  ADC  system. 
The  simulation  results  are  used  to  compare  the  performance  of  the  hardware  circuits  with 
the  ideal  performance.  Also,  the  design  and  analysis  of  a  parity  timing  processor  is  given. 
The  timing  processor  is  built,  and  the  clock  signal  of  the  processor  using  the  detected 
optical  pulse  signal  input  is  analyzed.  The  MATLAB  program  developed  in  this  thesis  can 
easily  be  applied  to  other  SNS  ADC  configurations  using  different  moduli 

One  of  the  goals  initially  envisioned  for  this  thesis  was  to  characterize  the  whole 
SNS  ADC  component.  This,  however,  had  to  be  changed  because  of  the  need  to  integrate 
the  folding  circuit  and  timing  processor.  Therefore,  this  thesis  discusses  only  the  static 
data  for  the  timing  processor  characteristics.  For  future  research,  signal  analysis  of  the  data 
bus  and  timing  processor  are  needed  to  observe  the  time  difference  between  both  signal 
paths.  Also,  timing  analysis  of  the  data  I/O  at  the  latching  states  using  the  optical  pulses  is 
needed. 
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APPENDIX  MATLAB™  SOURCE  CODE 


A.  Vth_calc.m 

%  Thesis  program  No.l  (Vth_calc.m) 

%  This  m-file  generates  normalized  threshold  intensity  /  voltage 
%  and  its  difference  between  the  neighbors.  In  this  program  a 
%  numbCT  of  moduli  X  is  needed  as  an  input. 

%  Written  by  Hiromichi  Yamakoshi 

clc,clg  %  input  the  number  of  modulus,  X 

X  =  inputCWhat  number  of  modulus  do  you  want?  Please  input  number :  ’); 
Drange  =X*(0: 1/X:  1 );  %  creates  input  voltage  index 

%  from  0  to  X 

ouq)ut  =  (sin((pi*Drange)/(2*X))).'^2;  %  apply  Eq.  3 

n  =  l:length(oulput)-l; 

diffD  =  diff(ouq)ut);  %  calculate  the  difference  between 

%  each  two  normalized  threshold 

plot(Drange,output,Drange,output,'o'),grid  %  start  plotting  process 
title(['Calculation  of  Normalized  Threshold  in  modulus  ',num2str(X)]) 
xlabeK'Input  Voltage  Index') 
ylabel(’Normalized  Detector  Ouq)ut’) 

diffD  =[0  diffD];  %  display  normalized  threshold 

disp(['Vth  of  mod',num2str(X),'  dVth  of  mod',niun2str(X)]) 
disp([output'  diffD*]) 

B .  thmtobin.m 

%  Thesis  program  No.2  (thmtobin.m)) 

%  Diis  m-file  generates  diermometer  code  and  converted  binary  code  with  input  of 
%  desired  number  of  modulus  X. 

%  Hiromichi  Yamakoshi 
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clc 

X  =  inputCWhat  number  of  modulus  do  you  want?  Please  input  number :  ’); 

X  =  X;  %  First,  using  X,  create  bit  information 

n=  1; 

while  X  >  2 

X  =  0.5  *  x; 

n  =  n+1;  %  n :  bit  information 

end 

dispCEWe  need  ',int2str(n),'  bits.']); 
disp('and') 

disp([’possible  comparator  information  of  moduli  ',int2str(X),'  is’]) 
comparator  =  zeros(X);  %  Next,  create  the  thermometer  code 

for  i=  1  :X  %  comparator :  thermometer  code  info. 

comparator(i,:)  =[ones(l,i-l),  zeros(l,X-i),i-l]; 
end 

compout  =  comparator(:,  1  :X- 1);  %  reshape  the  thermometer  code  output 

state  =  comparator(:,X); 

disp(rot90(comparator)) 

dispC  cf.  the  first  row  number  shows  comparator  states.') 
bin  =  zeros(X,n);  %finally  get  the  binary  code 

for  i  =  1:X;  %  to  create  binary  code,  we  uses  the  index 

m  =  1;  %  information  X 

p  =  i-l;  %  p  creates  index  starting  0  to  X-1 

if  p  =  0  %  create  bin.  infonnation  “0” 

bin(i,:)=bin(i,:); 

elre  %  create  bit  information  “1” 

whOe  p  >=2 

q  =  rem(p,2); 
bin(i,n+l-m)  =  q; 
p  =  fix(p/2); 
m  =  m+l; 
end 
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bin(i,n+l-m)  =  1; 
end 
end 

binO  =  flipb^in); 

binout  =  [rot90(state);iDt90(bin0)] ;  %  reshape  the  binary  code  output 

disp('also') 

dispCbinary  information  of  each  state  is’) 
disp(binout) 

dispC  cf.  the  first  row  number  shows  comparator  states.') 

C.  Lab_simu.m 

%  Thesis  program  No.3(Lab_simu.m) 

%  This  m-file  generates  the  number  of  folds  and  thermometer  code  of  each  sample 
%  point  of  each  modulus,  the  number  of  sample  point  between  two  thresholds  of 
%  each  comparator,  and  predicted  normalized  LabVBEW  ou^ut  of  each  san^le 
%  point. 

%  This  program  is  designed  for  8-bit  SNS  ADC  LabVIEW  circuit 

% 

%  Hiromichi  Yamakoshi 

%  First,  we  input  LAB  VIEW  setting  information. 

total_index  =  5000;  %  5000  samples  are  used  as  the  total  index. 

pt_cycle9  =  total_index/110;  %  1 10, 99,  and  90  are  used  as  numbCT  of  cycle  of 
pt_cyclel0  =  total_index/99;  %  each  modulus,  and  generates  the  number  of 
pt_cyclel  1  =  total_index/90;  %  sample  point  per  cycle  of  each  modulus. 

%  using  the  info  above,  generates  the  number  of 
pt_state  =  pt_cycle9/(2*9);  %  sample  point  between  two  thresholds  of  each 

%  comparator 
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%  now  input  the  number  of  sample  points  of  observation.  Now  we  choose  modulus 
%  9,10  and  1 1  for  8  bit  SNS  ADC.  We  use  256  dynamic  range  states  out  of  990. 

%  So  we  actually  need  the  index  of  sample  number,  0  to  646 

index=0:646; 

mod9folds  =  fix(index^t_cycle9);  %  generates  number  of  folds  for  each 
modlOfolds  =  fix(index4>t_cyclel0);  %  modulus  relates  to  the  sample  number  of 
modi  Ifolds  =  fix(index/pt_cyclel  1);  %  the  index  which  makes  8  bit  resolution 

%  next  we  number  the  ststes  of  each  modulus.  Each  modulus  has  the 
%number  which  index  fix>m  0  to  2*(number  of  modulus)  -1  in  one  cycle 
%  (Ex.)  they  have  18  states  ( 0  to  17)  in  mod  9  (mod  9  has  9  states) 

mod9states  =  fix((index-(mod9folds*pt_cycle9))/pt_state); 
modlOstates  =  fix((index-(modlOfolds*pt_cyclelO))^t_state); 
modllstates  =fix((index-(modllfolds*pt_cyclell))/pt_state); 

Drange  =  mod9folds*(2*9)  +  mod9states;  %  generates  dynamic  range 

%  Next,  we  create  the  thermometer  code-expression  of  each  modulus  (i,e, 

%  00 1 1 1 1 1 1 1  for  mod  9,  and  so  on) 

Ml=9;,  M2=10;,  M3=ll; 
codel=zeros(length(mod9states),Ml-l); 
code2=zeros(length(modl0states),M2-l); 
code3=zeros(length(modl  lstates),M3-l); 

%  formod9 

for  i=l  :length(mod9states); 
if  mod9states(i)  <  Ml 

codel(i,:)=[2eros(l,Ml-mod9states(i)-l)  ones(l,mod9states(i))]; 
else 

code  1  (i,:)=[zCTOs(  1  ,mod9states(i)-Ml )  ones(  1 ,2*Ml-mod9states(i)- 1 )]; 
end 
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end 

%  for  mod  10 

for  i=l  :length(modlOstates); 
if  modlOstates(i)  <  M2 

code2(i,:)=[zeros(l,M2-modl0states(i)-l)ones(l,modl0states(i))]; 

else 

code2(i,:)=[zeros(l,modl0states(i)-M2)cHies(l^*M2-modl0states(i)-l)]; 

end 

end 

%  for  mod  11 

for  i=l  :length(modl  1  states); 
if  modi  1  states(i)  <  M3 

code3(i,:)=[zeros(l,M3-modl  Istates(i)-l)  ones(l,modllstates(i))]; 
else 

code3(i,:)=[zeros(l^odllstates(i)-M3)  ones(l,2*M3-modl  lstates(i)-l)]; 
end 
end 

%  display  and  record  each  modulus'  thermometer  code  and  dynamic  range 
diary 

diq)('index_No.  dynamic_range  thermo_code  of  mod  9’) 
disp([index'  Drange'  codel]) 

disp('index_No.  dynamic_range  thermo_code  of  mod  10’) 
disp([index'  Drange'  code2]) 

di^('index_No.  dynamic_range  thermo_code  of  mod  11') 
disp([index'  Drange’  code3]) 

%  finally  we  get  ou^ut  voltage  related  to  the  index  No. 

mod9Vout  =  sin(pi*index/pt_cycle9).'^2'; 

modlOVout=sin(pi*index/pt_cyclelO).'^2’; 

modllVout  =  sin(pi*index/pt_cyclel  1).'^2'; 

dispCsimulated  output  signal  of  each  modulus  are  as  follows') 

dispC  index_No  dynamic_range  mod9  out  modlO  out  modi  1  out’) 

disp([index’  Drange'  mod9Vout  modlOVout  modi  1  Vout]) 
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diary  off 


D.  idlprty.m 

%  Thesis  program  No.4(i(ilprty.m) 

%  This  m-file  generates  the  threshold  value  of  each  parity  circuit  in  the  number  of 
%  modulus  with  the  least  significant  bit  (LSB)  The  number  of  modulus  in  LSB  and 
%  the  desired  percentage  of  decimation  is  required  to  get  the  threshold  voltage. 

%  This  program  is  designed  for  8-bit  SNS  ADC  LabVIEW  circuit 
% 

%  Hiromichi  Yamakoshi 

%  Thesis  project.  Jul  27  1994 

clc,clg  %  Input  the  required  information 

m=input(Tlease  input  the  number  of  the  least  significant  mod. :') 
n=inputCAnd  input  the  desired  percent  of  the  decimation  width :') 

Drange  =  m*(0:l/m:l); 

Vin  ^i*Drange/(2*m); 

Vlh  =  sin(pi*Drange/(2*m)).'^2;  %  First,  get  threshold  value  of  mod.m 

dVin=Vin(2)-Vin(l );  %  get  d(Vin)  that  is  used  to  get  %decimation 

dtmod=n*dVin^(X);  %  get  ±0.5*n(%)  of  d(Vin)  from  Vth 

tl=Vin-dtmod; 
t2=Vin+dtmod; 

tmod=[tl  ;t2] ;  %  create  parity  check 

tmod=ieshape(tmod,  1 ,2*lenglh(Vin));  %  input  value  and  that 

trnod=tinod(:,2:length(tttnod)- 1);  %  element  vector 

PVth=sin(tmod).''2; 

plot(tmod,PVth,tmodT*Vth,’x')  %  compute  and  graphically  display  the 
set(gca,’XTickLabels',[])  %  threshold  value  of  parity  circuit 

hold  on 

Vlhline=Vth'*ones(  1 4ength(  Vin)); 
tolanceline=PVth'*ones(l  ,length(Vin)); 
for  i=l:length(Vin) 

plot(Vin,Vthline(i,:),'c-.’,Vin,tolanceline,'r-') 


96 


end 

hold  off 
PVth' 


E.  sys_simu.m 

%  Thesis  program  No.5(sys_simu.m) 

%  This  m-file  generates  the  number  of  folds  and  diamometer  code  of  each  sample 
%  point  of  each  modulus,  the  number  of  sample  point  between  two  thresholds  of 
%  each  ccMnparattn:,  and  predicted  normalized  LabVIEW  ou^ut  of  each  san^le 
%  point.  In  addition  to  these,  this  m-file  also  returns  the  closest  threshold  level  and 
%  its  difference  for  each  sampling  point,  and  parity  bit  information.  The  desired 
%decimation  information  is  required  to  get  these  value  above. 

%  This  program  is  designed  for  8-bit  SNS  ADC  LabVIEW  circuit 
% 

%  Hiromichi  Yamakoshi 


%  enter  the  desired  %tolCTance 
n=input('Input  the  desired  %decimation  : '); 


%  First,  we  have  to  set  up  LAB  VIEW  setting  information. 


total_index  =  5000; 
pt_cycle9  =total_index/110; 
pt_cyclelO  =  total_index/99; 
pt_cyclell  =  total_index/90; 

pt_state  =  pt_cycle9/(2*9); 


%  5000  sample  s  is  used  as  die  total  index. 

%  1 10,99,and  90  are  used  as  number  of  cycle  of 
%  each  modulus, and  generates  the  number  of 
%  sample  point  per  cycle  for  each  moduli. 

%  using  the  info  above,  we 
%  get  the  #  of  point  per 
%  threshold  state  for  each  mod. 


% 


%  now  input  the  #  of  sample  points  of  observation.  Now  we  choose  modulus 
%  9,10,and  11  for  8-bit  SNS  ADC.  We  use  256  dynamic  range  states  out  of  990 
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%  So  we  actually  need  the  index  of  san:5)le  number,  0  to  646. 


% 

mdex=0:646; 

nQod9folds  =  fix(index^t_cycle9);  %  get  the  #of  folding 

modlOfolds  =  fix(index^t_cyclelO);  %  for  each  moduli 

modi  Ifolds  =  fix(index/pt_cyclel  1 );  %  related  to  index 

%  Next  we  number  the  ststes  of  each  modulus,  each  mod.  has  the 
%  numbCT  which  index#  from  0  to  2*(#of  mod)  -1  in  one  cycle 
%  (Ex.)  they  have  18  states  in  mod  9.  (mod  9  has  9  states) 

mod9states  =  fix((index-(mod9folds*pt_cycle9))/pt_state); 
modlOstates  =  fix((index-(modlOfolds*pt_cyclelO))/pt_state); 
modi  Istates  =  fix((index-(modl  lfolds*pt_cyclel  l))/pt_state); 

Drange  =  nKxi9folds*(2*9)  +  mod9states;  %  get  related  dynamic  range 

%  Next,  we  create  the  thermometer  code-expression  of 
%  each  modulus  (i,e,  (Xllllllllfor  mod  9,  and  so  on) 

Ml=9;,  M2=10;,  M3=ll; 
code  1  =zeros(length(mod9states),M  1  - 1 ); 
code2=zeros(length(mod  lOstates)  ,M2- 1 ); 
code3=zeros(length(modl  lstates),M3-l); 

%  for  mod  9 

for  i=l  :length(mod9states); 
if  mod9states(i)  <  Ml 

code  1  (i,:)=[zeiDs(  1  ,M1  -mod9states(i)- 1 )  ones(l  ,mod9states(i))] ; 
else  %if  mod9states(i)  >=  Ml 

codel(i,:)=[zeros(l,mod9states(i)-Ml)ones(l,2*Ml-mod9states(i)-l)]; 
end 
end 

%  for  mod  10 
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for  i=l  :length(modlOstates); 
if  modlOstates(i)  <  M2 

code2(i,:)=[zeros(l>I2-modl0states(i)-l)  ones(l,modlOstates(i))]; 
else 

code2(i,:)=[zeros(l,modl0states(i)-M2)ones(l,2*M2-modl0states(i)-l)]; 

end 

end 

%  for  mod  11 

for  i=l  :length(modl  1  states); 
if  modi  lstates(i)  <  M3 

code3(i,:)=[zeios(l,M3-modl  lstates(i)-l)  ones(l,modl  lstates(i))]; 
else 

code3(i,:)=[zeros(l,modllstates(i)-M3)ones(l,2*M3-modllstates(i)-l)]; 

end 

end 

%  display  and  record  each  modulus'  thermometer  code  and  dynamic  range 
diary 

disp('index_No.  dynamic  range  thermo_code  of  mod  9’) 
disp([index’  Drange’  codel]) 

disp('index_No.  dynamic_range  thermo_code  of  mod  10') 
disp([index'  Drange'  code2]) 

disp('index_No.  dynamic_range  thermo_code  of  mod  11') 
disp([index'  Drange'  code3]) 
diary  off 

%  Now  we  get  ou^ut  value  related  to  the  index  No. 

% 

mod9Vout  =  sin(pi*index/pt_cycle9).'^2; 
modlOVout=  sin(pi*index/pt_cyclel0).''2; 
modi  lVout=  sin(pi*index/pt_cyclel  1).'^2; 

diary 
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disp(['index_No.  Vin_  index  mod9out  modlOout  modi  lout*]) 
disp([mdex'  Drange'  mod9Vout'  modlOVout’  modllVout’]) 
diary  off 

%  then  we  get  the  difference  between  two  threshold 
%  and  threshold  value  of  each  mod.  (for  one  cycle) 

% 

tl=linspace(0,pi^*Ml  +l);,Vthl=sin(tl).''2;dVthl=diff(Vthl);  %  get  threshold  & 
t2=linspace(0,pi;2*M2  +l);,Vth2=sin(t2).'^2;dVth2=diff(Vth2);  %  its  diff.  for 
t3=linq)ace(0,pi,2*M3  +l);,Vth3=sin(t3).'^2;dVth3=diff(Vth3);  %  each  modulus, 
diary 

disp([’threshold  value  of  mod  9*]) 
disp(Vthr) 

disp([’threshold  value  of  mod  10']) 
disp(Vth2') 

disp([’threshold  value  of  mod  11']) 
disp(Vth3') 

disp(['the  diffCTence  between  two  threshold  of  mod  9']) 
dVthl* 

disp(['the  difference  between  two  threshold  of  mod  101) 
dVth2' 

disp(['the  difference  between  two  threshold  of  mod  111) 

dVth3' 

diary  off 

clear  tl  t2t3 

%  Begin  to  get  the  closest  threshold  value  as  for  each  sample  point 
%  after  the  following  process  we  get  the  closest  threshold  value  and 
%  its  difference  for  every  sample  point 

Adiff9Vout  =zeros(l,length(index));,Bdiff9Vout=zeros(l,length(index)); 
AdifflOVout=zeros(l,length(index));,BdifflOVout=zeros(l,length(index)); 
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Adiffl  lVout=zeros(Uength(index));3diffl  1  Vout=zeros(l,length(index)); 
diff9Vout  =2eros(l,length(index));,closeVth9  =  mod9states; 
difflOVout=zeK)s(l,length(index));,closeVthlO  =  modlOstates; 
diff  1 1  Vout=zeios(l,length(index));,closeVthl  1  =  modi  Istates; 

% 

%  A:  get  difference  between  Vout@index  pt.  &  just  below  Vth  (for  positive  slope), 
%  or  between  VouHSindex  pL  &  just  above  Vth  (for  negative  slope) 

%  B:  get  difference  between  VouttSindex  pt.  &  just  above  Vth  (for  positive  slope), 
%  or  between  Vout@index  pt.  &  just  below  Vth  (for  negative  slope) 

for  i=l:lengtfi(index) 

Adiff9Vout(i)  =mod9Vout(i)  -Vthl(nK)d9states(i)+l); 
Adiffl0Vout(i)=modl0Vout(i)-Vth2(modl0states(i)+l); 

Adiffl  1  Vout(i)=modl  1  Vout(i)-Vth3(modl  lstates(i)+l); 

Bdifi9Vout(i)  =mod9Vout(i)  -Vthl(mod9states(i)+2); 
Bdiffl0Vout(i)=modl0Vout(i)-Vth2(modl0states(i)+2); 

Bdiffl  1  Vout(i)=modl  1  Vout(i)-Vth3(modl  lstates(i)+2); 
end 
% 

%  Compare  A  &  B,  and  show  the  diff. value  from  which  Vth's  stage. 

% 

for  i=l  :length(index) 

if  abs(Adiff9Vout(i))  <  abs(Bdiff9Vout(i)) 
diff9Vout(i)  =  Adiff9Vout(i); 
else 

difif9Vout(i)  =  BdifE9Vout(i); 
closeVth9(i)  ^od9states(i)+l; 
end 

if  abs(AdifflOVout(i))  <  abs(BdifflOVout(i)) 
difflOVoutfi)  =  AdifflOVout(i); 
else 

difflOVout(i)  =  BdifflOVoutd); 
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closeVthKXi)  =modlOstates(i)+l; 


if  abs(Adiffl  1  Vout(i))  <  abs(Bdiffl  1  Vout(i)) 
diff  1 1  Vout(i)  =  Adiff  1 1  Vout(i); 
else 

difn  1  Vout(i)  =  Bdiffl  1  Vout(i); 
closeVthl  l(i)  =modl  lstates(i)+l; 
end 
end 

%  improve  the  effect  of  very  top/bottom  threshold  case 
%  (on  this  MATLAB  code,  I  treat  the  normalized  "O"/"!" 
%  as  a  threshold,  but  actually  there  are  not.) 

%  for  mod  9 

for  i=l  :length(mdex) 

if  mod9states(i)  =  0; 

diff9Vout(i)  =  Bdiff9Vout(i); 
closeVth9(i)  =  1; 
elseif  mod9states(i)=  9; 

diff9Vout(i)  =  Bdiff9Vout(i); 
closeVth9(i)  =  8; 
elseif  mod9states(i)  ==  8; 

diff9Vout(i)  =  Adiff9Vout(i); 
closeVth9(i)  =  8; 
elseif  mod9states(i)  =17; 

-  diff9Vout(i)  =  Adiff9Vout(i); 
closeVth9(i)  =  1; 
elseif  closeVth9(i)  =10; 

closeVth9(i)  =  8; 
elseif  closeVth9(i)  =11; 

closeVth9(i)  =  7; 
elseif  closeVth9(i)  =12; 
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closeVth9(i)  =  6; 
elseif  closeVth9(i)  =13; 

closeVth9(i)  =  5; 
elseif  closeVth9(i)  =14; 

closeVth9(i)  =  4; 
elseif  closeVth9(i)  =15; 

closeVth9(i)  =  3; 
elseif  closeVth9(i)  =16; 

closeVth9(i)  =  2; 
elseif  close Vth9(i)  =17; 

close Vth9(i)  =  1; 
end 
end 

%for  mod  10 
for  i=l  :length(index) 

if  modlOstates(i)  =  0; 

difflOVout(i)  =  BdifflOVout(i); 
closeVthlO(i)  =  1; 
elseif  modlOstates(i)=  10; 

difflOVout(i)  =  BdifflOVout(i); 
closeVthlO(i)  =  9; 
elseif  modlOstates(i)  =  9; 

difflOVout(i)  =  AdifflOVouKi); 
closeVthlO(i)  =  9; 
elseif  modlOstates(i)  =19; 

difflOVout(i)  =  AdifflOVout(i); 
closeVthlO(i)  =  1; 
elseif  closeVthlO(i)  =11; 

closeVthlO(i)  =  9; 
elseif  closeVthlO(i)  =12; 
closeVthlO(i)  =  8; 
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elseif  closeVthlO(i)  =13; 

closeVthlO(i)  =  7; 
elseif  closeVthlO(i)  =14; 

closeVthlO(i)  =  6; 
elseif  closeVthlO(i)  =15; 

closeVthlO(i)  =  5; 
elseif  closeVthlO(i)  =16; 

closeVthlO(i)  =  4; 
elseif  closeVthlO(i)  =17; 

closeVthlO(i)  =  3; 
elseif  closeVthlO(i)  =18; 

closeVthlO(i)  =  2; 
elseif  closeVthlO(i)  =19; 

closeVthlO(i)  =  1; 
end 
end 

%for  modi  1 

for  i=l  :length(index) 

if  modllstates(i)  =  0; 

diffl  1  Vout(i)  =  Bdiffl  1  Vout(i); 
closeVthll(i)  =  1; 
elseif  modi  lstates(i)=  11; 

diffl  1  Vout(i)  =  Bdiffl  lVout(i); 
closeVthll(i)  =  10; 
elseif  modi  1  states(i)  =10; 

diffl  lVout(i)  =  AdiffllVout(i); 
closeVthll(i)  =  10; 
elseif  modllstates(i)  =21; 

diffl  lVout(i)  =  Adiffl  lVout(i); 
closeVthll(i)  =  1; 
elseif  closeVthll(i)  =12; 
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closeVthll(i)  =  10; 
elseifcloseVthll(i)==13; 

closeVthll(i)  =  9; 
elseifcloseVthll(i)=14; 

closeVthll(i)  =  8; 
elseif  closeVthl  l(i)  =15; 

closeVthll(i)  =  7; 
elseif  CloseVthl  l(i)  =16; 

closeVthll(i)  =  6; 
elseif  CloseVthl  l(i)  =17; 

CloseVthl  l(i)  =  5; 
elseif  closeVthl  l(i)  =18; 

closeVthl  l(i)  =  4; 
elseif  closeVthl  l(i)  =19; 

closeVthll(i)  =  3; 
elseif  CloseVthl  l(i)  =20; 

closeVthll(i)  =  2; 
elseif  CloseVthl  l(i)  =21; 

closeVthll(i)  =  1; 
end 
end 
diaiy 

disp([index  Drange'  diff9Vout'  closeVth9'  difflOVout'  closeVthlO'  diffllVout' 
CloseVthl  n) 
diary  off 

clear  Adiff9Vout  AdifflOVout  Adiffl  1  Vout  Bdiff9Vout  BdifflOVout  Bdiffl  1  Vout 

%  and  finally,  create  parity  circuit  portion 

t  =  linspace(0,pi/2,Ml+l); 

dt=t(2)-t(l); 

dtmod=n*dl/200 

tlow=t-dtinod;,  tup=t+dtmod; 
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tm(xi=[tIow;tup];,tmod^eshape(tmod,l,2*length(t));  %create  the  threshold 

tmod=tmod(:^:length(tmod)-l);  %  value  of  pari 

Vparity=sin(tmod).'^2; 

Vpaiity=reshape(Vparity,2,Ml); 

Vprtyup=Vparity(  1  ,:);,Vprtydn=Vparity(2,:); 

Vprtypitpdiff=Vprtyup-Vthl  (:,1 :9);,Vprtypupdiff=Vprtypupdiff(2:9); 
Vprtypdndiff=Vprtydn-Vthl(:,l:9);,Vprtypdndiff=Vprtypdndiff(2:9); 
Vprtyniqxliff=Vlhl  (:,2: 10)- Vprtyup;,Vprtynupdiff=Vprtynupdiff(  1:8); 
Vprtyndndiff=Vthl  (:,2: 10)- Vprtydn;,Vprtyndndiff=Vprtyndndiff(  1:8); 
pityind=zeros(  1 4ength(index)); 

for  i=l  :length(index);  %create  parity  bit  infomaation 

ifdiff9Vout(i)<0 

if  abs(diff9Vout(i))  >  abs(Vprtyndndiff(closeVth9(i)))  & ... 
abs(diff9Vout(i))  <  abs(Vprtynupdiff(closeVth9(i))) 
prtyind(i)=0; 
else 

prtyind(i)=l; 
end 


%  value  of  parity 


if  abs(diff9Vout(i))  >  abs(Vprtypupdiff(closeVth9(i)))  & ... 
abs(diff9Vout(i))  <  abs(Vprtypdndiff(closeVth9(i))) 
prtyind(i)=0; 
else 

prtyind(i)=l; 

end 


diary 

disp([index'  Drange'  diff9Vout’  closeVth9’  p^tyind^) 
diary  off 
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